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ABSTRACT 


A  theoretical  and  experimental  investigation  is 
presented  for  conduction  through  and  free  convection 
from  a  tapered,  downward-pro jecting  fin  immersed  in  an 
isothermal,  quiescent  fluid.  The  problem  is  solved  by 
assuming  quasi-one-dimensional  heat  conduction  in  the 
fin  and  matching  the  solution  to  that  of  the  convection 
system,  which  is  treated  as  a  boundary  layer  problem. 

The  theoretical  analysis  undertaken  for  an  infinite 
Prandtl  number  reveals  a  solution  which  takes  the  form 
of  a  power  law  temperature  distribution  along  the  fin. 

It  is  shown  that  this  power  (n)  is  independent  of  the 
precise  form  of  the  fin  profile  and  is  dependent  only 
on  a  single  non-dimensional  group  The  effect  of  n 

on  heat  transfer,  skin  friction,  fin  effectiveness  and 
the  corresponding  boundary  layer  profiles  is  discussed. 

In  comparing  these  results  to  those  from  previous  work 
for  the  unit  Prandtl  number  case,  it  is  demonstrated 
that  the  independent  effect  of  Prandtl  number  is  very 
small  for  the  range  1  <1  a  <L  co  . 

Theoretical  predictions  were  very  well  verified  by 
the  experimental  results  which  consist  of  temperature  pro¬ 
files  and  the  n  -  x  relation  for  different  fin  geometries 
and  surrounding  fluids. 
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NOMENCLATURE 

X 

vertical  distance  above  fin  tip 

y 

lateral  distance  from  fin  surface 

w 

fin  length 

6 
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n 

@ 
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dissipation  function 

1 
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ix 


H,T,Q,R  functions  defined  in  text 

a,b,b' ,c,A,d,g  coefficients  in  power  series 


a, 7 

eigenvalues 

r 

gamma  function 

k 

thermal  conductivity 

h 

heat  transfer  coefficient 
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non-dimensional  group  defined  in  text 
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fin  effectiveness 
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body  force  per  unit  mass 
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f 
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CHAPTER  I 


INTRODUCTION 


It  is  well  known  that  extending  the  basic  heat  trans¬ 
fer  surface  by  means  of  fins  considerably  increases  the 
rate  of  heat  transfer  from  a  surface  immersed  in  an  iso¬ 
thermal,  quiescent  fluid.  The  quantity  most  frequently 

used  in  evaluating  the  performance  of  fins  is  fin  effect- 

* 

iveness,  In  the  past  1  ,  fin  effectiveness  has  been 

— 

based  on  the  solution  of  the  conduction  problem  with  the 
convective  boundary  condition  given  by  an  empirical,  con¬ 
stant  value  of  heat  transfer  coefficient. 

Generally  the  surface  heat  transfer  coefficient  is 
not  constant  but  is  some  function  of  distance  along  the 
fin  surface  and  therefore  particular  attention  must  be 
given  to  its  form.  For  a  vertical  fin  as  considered  in 
the  following  analysis,  the  problem  may  be  solved  by 
treating  it  as  a  quasi-one-dimensional  conduction  problem 
with  the  convective  boundary  condition  obtained  from  the 
solution  of  the  boundary  layer  problem.  First  this  en¬ 
tails  solving  the  problem  of  laminar  free  convection  from 
a  vertical,  nonisothermal  surface.  Secondly,  the  surface 
temperature  and  heat  flux  in  the  corresponding  conduction 


Numbers  designate  references  given  on  page  57  . 


' 


problem  is  matched  to  that  obtained  in  the  convection 
problem. 
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In  most  practical  applications  the  surface  from 
which  heat  transfer  occurs  is  nonisothermal .  There  is 
a  very  limited  amount  of  information  available  in  the 
literature  for  free  convection  from  a  vertical,  noniso= 


thermal  surface.  Sparrow  and  Gregg 


give  results  for 


the  case  of  uniform  heat  flux,  based  on  numerical  solu- 
tions  of  the  differential  equations  of  the  boundary  layer 


The  above  authors 


also  considered  the  power  law  and 


exponential  forms  of  surface  temperature  for  the  case  of 


gases,  i.e.,  1. 

In  the  investigation  to  follow,  the  problem  of  con¬ 
duction  through  and  laminar  free  convection  from  a  tapered, 
downward-projecting  fin  with  a  large  slenderness  ratio, 
w/6r,  illustrated  in  figure  1.1,  is  considered.  A  large 
slenderness  ratio  allows  one  to  neglect  the  effects  of 
curvature  and  inclination  when  solving  the  flow  system 
as  a  boundary  layer  problem.  Both  theoretical  and  experi¬ 
mental  results  for  a  single  fin  are  presented;  in  particular 
for  the  case  in  which  the  fin  surface  temperature  varies 
in  a  simple  power  law  form.  The  effects  of  slenderness 
ratio  and  surrounding  fluid  were  investigated.  For  the 
range  of  conditions  considered  there  was  good  overall  agree¬ 
ment  between  theory  and  experiment. 


. 


3 


6  =  6r  (x/w)m 


FIG.  1.1  Tapered,  Downward-Projecting  Fin 


■ 
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Several  experimental  investigations 


4,5, 6, 7 


of 


free  convection  from  fins  have  been  undertaken  in  recent 
years.  Unfortunately,  these  were  primarily  concerned 
with  the  overall  behavior  of  banks  of  fins  under  varying 
conditions  and  their  results  cannot  be  directly  compared 
with  those  obtained  in  this  thesis. 

The  overall  problem  as  posed  in  this  thesis,  is  that 
of  a  slender,  downward-pro jecting ,  vertical  fin  of  known 
material,  shape  and  root  temperature  immersed  in  a  given 
isothermal,  quiescent  fluid. 
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CHAPTER  II 


THEORETICAL  ANALYSIS 


2.1  BOUNDARY  LAYER  PROBLEM 


An  analysis  is  made  for  laminar  free  convection  from 
a  downward-projecting  fin  having  a  non-uniform  surface 
temperature . 

2.1-1  FORMULATION  OF  PROBLEM 


The  equations  expressing  the  conservation  of  mass, 

* 

momentum  and  energy  for  laminar  incompressible  flows  due 
to  heat  conduction  from  a  surface  to  a  fluid  in  a  uniform 

=  0 

—  1  9  ^ 

~  VP  +  VV  V  +  G 

P 

K  v2t  + 

pcp 

Consider  the  application  of  these  equations  to  two 
dimensional  steady  flow  in  the  semi-infinite  region  illus¬ 
trated  in  figure  1.1.  On  applying  these  equations,  certain 
simplifying  approximations  are  made  with  respect  to  fluid 
properties.  All  properties  except  density  may  be  taken 


body  force  field  ares 


MASS 


V  •  V 


MOMENTUM: 


D\f 

Dt 


ENERGY: 


DT 

Dt 


* 


Following  Boussinesq. 


' 


6 


as  constant  when  dealing  with  moderate  temperature  changes 


8 


The  density  is  considered  as  a  variable  function  of 


temperature  only  in  forming  the  buoyancy  term.  The 

Boussinesqian  approximations  i.eOJ  — c^.  1  and  p  -  p 

^oo 

are  invoked  along  with  the  coefficient  of  expansion,  f3, 
defined  from  the  equation  of  state. 

Since  it  is  known  a  priori  that  a  boundary  layer  is 
formed  there  will  be  a  resulting  change  in  pressure  and 
the  sign  of  this  departure  from  the  hydrostatic  pressure 
depends  upon  whether  the  fin  is  being  heated  or  cooled. 
Taking  the  pressure  as  P  =  Pg  +  P^,  equations  (1)  may  be 
written  as: 


oo 


f  0 


d  u  ,  d  u 

u  -v - r  v  -v— 

ox  oy 


du  .  dv  _  n 
Sx  +  'Jy  ~  0 


1 

P 


^d 


d  v  ,  d  v 

u  t - r  v  3— 

ox  oy 


+  V 


1  dpd 

p  ~Sy 


'd2u  +  afu\  +  pg(T 
1  dy  ' 

+  v  flfv  +  aV 


T  ) 
00 


ix 


dy 


(2) 


dT  ,  dT 
u  +  v 
ox  dy 


=  K. 


d^T  d^T' 


+ 


2v 


ix 


dy 


'du 

.35 


+ 


IV 


2-1 


and  are  subject  to  the  following  boundary  conditions: 

y  =  0  (x  >  0)  u  =  v  =  0  Tq  =  Tq(x) 

y  =  00  u  =  v  =  0  T  =  T^ 


(3) 


. 
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In  accordance  with  the  usual  boundary  layer  approxi- 
mations,  y<<x,  and  hence  from  the  continuity 

equation  v«u„  Using  these  approximations  and  neglecting 
dissipation  a  posteriori ,  equations  (2)  become 


du  ,  dv  _ 

^  +  37  _  0 


>u 


u  S5E 


du  _ 

v5?  - 

dT  , 

U  37  +  V  3y 


Pg(T 

5t 


Toq)  +  V 


.2 

o  u 

^y2 


(4) 


/c 


b2T 

dy2 


Equations  (4)  may  be  written  in  terms  of  a  stream 
function s  f  =  tf/(x9y)  such  that  u  =  ^  and  v  =  -  thus 

satisfying  the  continuity  equationc  The  velocity  components 
u  and  v  are  thus  replaced  by  the  derivatives  of  the  stream 
function  and  from  this  substitution  there  results  a  pair 

s 

of  partial  differential  equations  for  ^  and  T  as  functions 
of  x  and  ye 

dK,  a2*  d*  aj^  _  pg(T  _  T  ^ 


oil/  a  ih  dib 

3y  dxoy  §x  ^  2 


bti  dT 


T  )  +  V  0 

00  ay3 


(5) 


37 


aT  _ 

57  ' 


/c 


d2T 

dy2 


2  0 1-2  SIMILARITY  SOLUTIONS 


It  is  desirable  to  seek  a  similarity  transformation 


which  will  transform  the  partial  differential  equations  to 


*  o,  =  £a*«i. 

p 

**  For  a  complete  analysis  of  reduction  of  equations  see 


' 


8 


a  pair  of  ordinary  differential  equations.  Sparrow  and 


Gregg 


carried  out  the  first  similarity  analysis  for 


f ree-convection  from  non-isothermal  surfaces.  Yang 
later  did  a  similar  analysis  and  in  both  papers  the  power 
law  and  exponential  temperature  distributions  were  dis¬ 


cussed.  It  appears 


10 


that  the  latter  is  not  really  a 


distinct  form  and  therefore  an  analysis  in  which  the  sur¬ 
face  temperature  takes  the  form 


0 


9  (” 

R  V  w 


n 


(n  >  0)  ,  (9  =  T  -  T  ) 

—  ’  co 


will  describe  the  general  similarity  situation  for  a  fin. 
This  situation  is  also  consistent  with  a  finite  tempera¬ 
ture  at  the  fin  tip  equal  to  that  of  the  surrounding  fluid 
Therefore,  the  similarity  variables  for  this  case  are 


taken 


as 


^(x,y) 


4v 


p£g  0R 

,2  n 
L4v  w  - 


rrf  3 


x 


f  (t)) 


©(x,y)  -  0 


R  (w 


n 


(rj) 


where 


■n 


rPg  © 


R 


a  2  n 
!-4v  w  - 


h 


1-n 


x 


These  variables  describe  the  general  case  but  the 
resulting  ordinary  differential  equations  contain  the 
Prandtl  number.  A  new  similarity  transformation  may  be 


(6) 


(^T  -  T  ■  0) 


;:r  3  >r  ua  ©  tt  to  rffiii  .  To  *  q  *  r: is  orfi  uzui 


found  which  will  render  the  equations  independent  of 
Prandtl  number.  This  is  done  by  considering  the  asymp¬ 
totic  situations  of  a  — and  oo .  For  this  analysis  the 


9 


situation  of  a — >-oo  will  be  dealt  with  and  will  give  re¬ 
sults  which,  as  well  as  representing  an  upper  bound,  can 
be  expected  to  lie  sufficiently  close  to  the  unit  Prandtl 
number  results  to  permit  accurate  interpolation. 

A  similar  transformation  to  that  utilized  by 


Lef evre 


11 


for  an  isothermal  surface  enables  the 


variables  for  a  non-isothermal  surface  to  be  written  as 

n+3 

x  4  F(£) 


^(x,y)  =  4/c 


Pg  e, 


L4/c2wn(l+a)  - 


€>(x,y) 


x 


n 


tw ) 


where 


i  = 


(3g  © 


R 


-4/c2  wn  ( 1+a)  - 


h 


1-n 


x 


for  the  extreme  values  of  Prandtl  number 


2.1-3  GOVERNING  EQUATIONS 

To  obtain  a  set  of  ordinary  differential  equations, 
the  similarity  variables  (7)  are  substituted  in  the 
boundary  layer  equations  (5)  to  give 


(7) 


- 


10 


I  I  0 


+  $ 


0  0 


+  (n  +  3)  FF  -2(n+l)(F)  +<J>  =  0 


i  i 


(8) 


0  +  (n  +  3)  FO 


4n  F  0  =  0 


These  equations  become  independent  of  Prandtl  number 
for  the  extreme  values  of  a  =  0  and  oo .  In  the  latter  case 

i  i  i 

the  momentum  equation  reduces  to  F  =  and  equations  (8) 


reduce  to 


V  T\7  '  1,1 

F  +  (n  +  3)  FF  -  4n  F  F  =  0 


(9) 


in  combined  form.  The  associated  boundary  conditions  ares 


•  i  i 


F  (o)  =  F  (o)  =  F  (o)  +  1  =  0 

ii  iii 

F  (oo )  =  F  (oo )  =  0 


(10) 


These  differ  from  the  unit  Prandtl  number  boundary 

i  • 

conditions  in  just  one  respect;  that  is,  F  (oo)  =  0 
■ 

replaces  F  (oo)  =  0  which  is  regarded  as  less  admissible 


as  o — »-oo  . 


2.2  SOLUTIONS  TO  THE  BOUNDARY  LAYER  EQUATIONS 


The  combined  equation  (9)  may  be  solved  by  the  method 


of  Meksyn 


12,13 


which  will  give  the  heat  transfer  and 


skin  friction  results.  From  these  results  the  velocity 
and  temperature  distributions  may  then  be  determined  by 
a  forward  integration  technique. 


o  -  1  +  <°J  ’  -  -  <°>*  .  ij 

I  I 


1  fsi  \  tollman.  sno  ^ajjfc  flx  anoirfionoo 


o  i 


11 


2.2-1  MEKSYN’S  METHOD 


In  obtaining  solutions  to  equation  (9) ,  use  is  made 

1  1  1  XV 

of  the  fact  that  F  and  F  are  large  only  in  a  narrow 
region  close  to  the  fin  surface.  This  fact  makes  it 
possible  to  seek  asymptotic  solutions  of  the  equations. 

Formally  solving  equation  (9)  as  a  linear  non- 
homogeneous  first  order  equation  in  F^"V  gives 

F 


IV 


-HU) 


FIV(o) 


4 


4n  F  F 


a  l  I 


HU) 


ai- 


ui) 


where 


H(£)  =  (n  +  3)  F  df* 


A  series  solution 


ao  *r 

a? 


,  F(?)  =  £ 


is  now  assumed 


r  =  o 


which  allows  (11)  to  be  written  in  the  form 


where 


r  =  o 

is  regarded  as  a  slowly  varying  function  of  £• 

Expanding  about  the  origin  in  a  Maclaurin  series  and 
utilizing  the  differential  equation  and  boundary  conditions 
there,  the  a^  are  readily  found  in  terms  of  the  parameter 

•  «  jy 

n  and  the  eigenvalues  F  (o)  and  F  (o) ,  i.e.,  a  and  y9  say. 
Expansion  of  (12)  and  comparison  of  terms  in  equal  powers  of 


(12) 

(13) 


r  io  qa  ni  (  -  .  no  aoaqxS 


12 


and 


the  b 

r 

in  terms  of 

the  a 

r 

and  consequently  in 

and  7 

« 

The  first 

few  terms  of  a  and  b  ares 

r  r 

a 

o 

= 

0 

a5 

=  0 

ai 

ss 

0 

a6 

=  -4n  a 

a2 

= 

a 

a7 

-  4n  +  (7n  -  3)  7a 

a3 

= 

i 

• 

O 

a8 

=  (3  -  15n)  7 

a4 

— 

7 

b 

o 

= 

7 

bl 

= 

0 

b2 

= 

-4n  a 

b3 

— 

4n(l  +  2ay) 

b4 

= 

~16n  7 

b5 

= 

16n  (72  - 

_  2 

5na  - 

3a2) 

b6 

= 

316n2a  + 

220n2a27 

2 

+  180na  7  +  180na 

To  facilitate  the  integration  of  equation  (12) ,  £  is 

replaced  by  a  new  independent  variable  t  defined  as 

co 


T  =  H(5)  -  i 


c  i' 

r  ^ 


r  =o 


where  the  c^  are  found  from  the  expansion  of  H(|)  giving 


the  first  few  terms  ass 

(n  +  3)  a, 


3! 

(n  +  3)  a. 


(n  +  3)  a( 
— 

(n  +  3)  a. 
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8: 


■ 

13 


(n  +  3) 


(n  +  3)  a 


8 


9! 


c3  =  0 


It  is  now t necessary  to  express  £  in  terms  of  t  by 
inverting  the  series  for  t,  taking 


« ■  i 

s  =  o 


00  A  t^1)/3 
s 


s  +  1 


(14) 


which  is  valid  for  sufficiently  small  values  of  t.  Follow¬ 


ing 


12 


,  Ag  is  the  coefficient  of  |  in  the  expression 

<co  +  cx4  +  c2i2  +  c343  +  ..J-(s+1)/3 


giving 


An  =  c 
0  o 


-1/3 


A. 


-  3  C1  Co 


5/3 


A. 


A3 


=  c 


-4/3 


28 

9 


C1C2 


140 

81 


'Cl\3n 


5/3 


5  C4  ,20 

1  ~  +  ~ 
o 


C2\2 


2 

220  C1  C2  ,  770 

27  3  243 


'C  \  4' 
C1V 


etc . 

The  A  9  through  the  c  and  hence  the  a  ,  implicitly  involve 
n,  a  and  y. 

Formally  integrating  equation  (12)  once  and  using  the 

i  i  i 

boundary  condition  at  the  origin  on  F  gives 


0' 


5  3o  sin  r.5  J  >oo  IJ  i  £ 


- 


14 


i  i  i 


e_T  *-<?>  f; 


dT  -  1 


(15) 


Utilizing  equations  (13)  and  (14)  and  re-writing  R(£) 


di 

dx 


as 


R(?) 


d£  =  t“2/3 

dT  T 


oo 


i 


s  =o 


d  t 
s 


s/3 


(16) 


where  3d  is  the  coefficient  of  £  in 
s 


+  +  2  l  +  ~1T  +  *  •  *XCo  +  Cl^  +  c2^ 2  +  c3^ 3  +  •  • 


3! 


(s+1) 


gives 


1  c  -1/3  b 
3  o  Do 


4  c  ~2/3 

3  o 


-  ~  b  - 
3  o  c 


Cln 


o- 


1  -1 

*=r  c 
3  o 


rb. 


bici 


c,\2n 


-  b  -2  +  b  P") 

o  c  o \c  / 

o  \  o7 


etc . 


Substituting  (16)  in  (15) ,  we  obtain 


i  i  i 


oo 


s  =o 


r  -t  T  ( s+i)  /3  - 1  dT  _ , 


Jo 


I  I  I 


so  that  the  condition  F  (oo)  =  0  is  satisfied  if 

oo 


I 


a. r  m  ■  1 


(17) 


s  =o 


*  For  brevity  and  due  to  their  complexity  the  remaining 

coefficients,  d  ,  are  omitted  here.  They  are  given  in 

s 

the  computer  program  in  Appendix  A,  p.  59  „ 


(?■ ) 
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This  series  gives  only  one  equation  for  the  eigen¬ 
values  a  and  7  in  terms  of  the  parameter  n  and  therefore 
a  second  equation  is  required.  This  second  equation  is 

1  1 

obtained  by  satisfying  the  condition  F  (00)  =  0  in  a 

'  '  1  -H ( s ) 

similar  manner  by  letting  F  (£)  =  e  s  Q(|)  ana  pro¬ 
ceeding  as  above.  Q(*)  results  in  a  new  set  of  values, 

1 

b  r,  the  first  few  being 


b 


o 


b 

b 

b 

b 

b 


1 

2 

3 

4 

5 


b 


6 


-1 

7 

0 

-a(5n  +  3) 

5n  +  3  +  ayflln  +  9) 

-7 (2 In  +  15) 

21n72  +  1572  -  130n2a2  -  228na2  -  90az 


The  a  .  c  and  A  terms  are  as  oreviously  obtained 
r  ’  r  s 

and  the  d,  are  replaced  by  g„  which  incorporate  the  new 
s  s 

■ 

values  b  This  results  in 


CD 


S  =0 


(18) 


where  the  g  ,  like  the  d  ,  contain  a,  7  and  the  temperature 
J  s  *  s 

index  n.  Therefore  equations  (17)  and  (18)  contain  the 
two  unknown  eigenvalues  and  constitute  the  essential  solu¬ 
tion  of  the  problem. 


. 
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2.2-2  SOLUTION  FOR  EIGENVALUES 

Since  the  simultaneous  equations  (17)  and  (18)  contain 
a  and  7  implicitly  it  is  convenient  to  solve  them  numeri¬ 
cally  by  an  iterative  procedure,  following  truncation  after 
seven  terms.  This  procedure  was  carried  out  by  a  computer 
program  as  given  in  Appendix  A,  written  in  Fortran  IV  langu¬ 
age.  The  calculations  were  carried  out  on  an  IBM  7040-1401 
machine  in  the  Department  of  Computing  Science,  University 
of  Alberta. 

The  basic  procedure  involved  in  the  program  is  given 
as  follows.  The  program  was  initially  set  up  for  seven 
terms  in  each  of  the  series.  With  the  temperature  index 
n  =  0,  an  initial  value  of  both  a  and  7  was  assumed.  Each 
series  was  then  summed  and  checked  with  the  correct  sum, 
i.e.,  unity  for  equation  (17)  and  -a  for  equation  (18). 

If  these  sums  were  not  within  the  set  limit  of  their  correct 
value,  a  was  incremented  by  an  amount  proportional  to  the 
error  in  one  series  while  7  was  held  constant.  Using  the 
new  value  of  a  the  series  were  again  summed  and  if  found 
to  be  still  in  error,  7  was  incremented  by  an  amount  pro¬ 
portional  to  the  error  in  the  second  series  while  a  was  held 
constant.  This  iterative  procedure  was  repeated  until  the 
two  series  came  within  the  set  limit  of  their  correct  sum 
thus  producing  acceptable  values  of  ct  and  7  for  n  =  0.  The 


' 
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procedure  was  repeated  for  the  range  of  values  of  n  from 
0  to  3  in  increments  of  0.1.  In  each  case  the  initial 
estimates  of  a  and  y  used  to  begin  the  iterative  procedure 
were  those  accepted  for  the  previous  value  of  nQ 

In  summing  each  of  the  series  (17)  and  (18)  frequent 
use  was  made  of  Euler's  transformation  to  remove  divergence 
when  necessary.  Since  many  summations  were  required  in  the 
iterative  process  it  was  necessary  to  adopt  a  systemmatic 
application  of  the  Euler  transformation:  this  was  incor¬ 
porated  in  a  subroutine  to  the  main  program. 

When  the  series  exhibited  divergence,  the  transforma¬ 
tion  was  begun  at  the  first  term  but  additional  transforma¬ 
tions,  each  starting  with  a  progressively  higher  term  were 
also  carried  out.  If  none  of  the  resulting  series  were 
convergent  the  process  was  repeated  on  the  first  transforma 
tion.  However,  if  one  or  more  of  the  first  set  of  transfor 
mations  gave  convergence  the  series  with  the  smallest  last 
term  was  used  in  the  analysis,  which  then  proceeded  in  the 
iterative  process.  In  the  majority  of  cases  the  series 
chosen  in  this  manner  also  resulted  in  the  series  with  the 
largest  number  of  unchanged  terms. 

The  eigenvalues  a  and  y  thus  obtained  are  plotted 
versus  n  in  figure  2.1  together  with  the  corresponding 


values  of  y  for  a  =  1 


tained  by  other  authors. 


and  the  extreme  values  of  y  ob* 


* 

fUfi/n  erf^  oi  sniijjo^dxra  b  ni  boiffi'xoq 


:i.  > 

.atoiilue  uorfio  vcf  banijs:* 


Heat  Transfer,  7  or  Skin  Friction 


18 


Temperature  Exponent  n 


Fig„  201  Heat  Transfer  and  Skin  Friction  Results 
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2.2-3  VELOCITY  AND  TEMPERATURE  PROFILES 

Following  Meksyn,  the  velocity  and  temperature  pro¬ 
files  can  be  expressed  in  terms  of  incomplete  gamma  func¬ 
tions  and  their  integrals  but  it  is  much  more  convenient 
here  to  use  a  and  y  in  a  forward  integration  procedure „ 

The  profiles  were  obtained  by  another  computer  pro¬ 
gram  also  written  in  Fortran  IV  and  given  in  Appendix  B. 
This  program  simply  employed  a  normal  forward  integration 
procedure,  in  this  case  using  a  fourth  degree  polynomial 

.  '  • • •  iv 

with  increments  of  0.05  to  determine  F,  F  ,  F  and  F 
V 

with  F  being  determined  from  the  differential  equation. 

ii  iii 

Since  the  conditions  on  both  F  and  F  at  infinity  are 
satisfied  in  Meksyn ' s  procedure,  it  was  only  necessary 

•  i  • 

to  forward  integrate  until  F  approached  zero.  This 
also  provided  a  further  check  on  the  eigenvalues  since 

i  i 

F  should  approach  zero  at  approximately  the  same  value 

i  i  i 

of  £  as  does  F  .  The  profiles  are  shown  in  figure  2.2 
which  illustrates  the  effect  of  varying  the  temperature 
index. 

2.3  MATCHING  CONDUCTION  PROBLEM 

Heat  transfer  within  the  fin  is  considered  as  a  quasi- 
one-dimensional  conduction  problem  with  the  heat  transfer 
coefficient  and  temperature  at  the  fin  surface  varying  in 


1  ;  >*.  cCoicj  nox3oi/Jbno3  t  10  e n©..?j  j ~-onc 
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FIGo  2  o  2  Velocity  and  Temperature  Profiles 
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the  manner  prescribed  from  the  solution  of  the  free  con¬ 
vection  problem  discussed  above 0 


2 „ 3-1  SOLUTION  OF  CONDUCTION  EQUATION 


The  equation  describing  quasi-one-dimensional  con¬ 
duction  within  the  fin  is 


d 

dx 


h  9 

o 


(19) 


where  the  fin  profile  6(x)  has  been  prescribed  and  the 
form  of  the  surface  heat  transfer  coefficient  h(x)  is  known 
from  the  boundary  layer  solution 0  Thus 


h  = 


(i 


|0 

■f  v^A.=o 

15 


i  „  e 


h  = 


kf  7 


r  f3g0 


R 


n-1 


■4/c^  wn(l+a)- 


x 


(20) 


Substituting  (20)  along  with  the  fin  profile,  6  =  6R 
into  the  conduction  equation  (19)  gives 

n-1 


I 


.m 


d20  ,  d9 

m  o  ,  m-1  o 

x  - x-  +  mx 


dx4 


dx 


w 

^sA5R. 


m\  r 


Pg  ©R  -\h 


L4/c2wn(l+a)  J 


7  x  0=0 


(21) 


The  fin  solution  is  compatible  with  the  boundary  layer 
solution,  only  if  (21)  has  the  solution  0-  =  0R 

In  the  conduction  problem  any  power  law  temperature 
solution  requires  that  the  exponent  of  x  in  the  coefficient 
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of  the  second  derivative  is  two  greater  than  that  of  the 

corresponding  coefficient  of  the  function  0  .  That  is. 

o  ’ 

n  =  4m  -  7  .  (22) 


The  indicial  equation  gives  another  relation  between 
n  and  m;  i.e. 

r 


n  =  -|(1  -  m)  +  -| 


'■  - 1,2  * 


/_  m\r  Pg  9 


R 


L4/c2wn  (1+cr)  -1 


k 


y 


(23) 


where  the  positive  root  is  taken  to  satisfy  the  requirement 
that  n  >  0.  The  elimination  of  m  from  (22)  and  (23)  results 

H 


m 


n 


2  +  J_ 

10  10 


where 


X  "  1 7“ 


w 

AAV 


.  ,  80 
9  +  —  X7 
\f2 

-Pg  eR  w3-i^ 
-/c2  (1+a)  - 


(24) 


(25) 


Thus  there  is  only  one  value  of  the  dimensionless 
group  x  f°r  value  of  n,  for  which  the  power  law  fin 

surface  temperature  profile  is  valid  within  the  approxi¬ 
mations  of  this  analysis.  The  relation  between  n  and  x  f°r 
a  =  oo  is  illustrated  in  figure  2.3  along  with  corresponding 
theoretical  points  for  a  =  1. 


2.3-2  FIN  EFFECTIVENESS 


Although  equation  (24)  gives  the  solution  to  the  fin 
conduction  problem  in  terms  of  temperature  index  n,  it 
does  not  fully  describe  the  heat  transfer  characteristics 
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of  the  fin  itself .  The  quantity  that  best  describes  the 
fin  behavior  from  the  heat  transfer  aspect  is  the  fin 


effectiveness,  normally 


defined  as 


>w 

@q(x)  dx 


eRw 


This  definition  however  assumes  that  an  empirical, 
constant  value  heat  transfer  coefficient  is  available 
but  as  this  analysis  illustrates,  h  is  not  constant  but 
is  a  variable  function  of  x0  Therefore,  the  above  de¬ 
finition  is  inadequate  for  present  purposes  and  will  be 

replaced  by  TW  dx 

€  =  "2 _ ° - 

hRSR  w 

where  the  “"reference"  heat  transfer  rate  is  taken  at  root 
conditions .  Substituting  for  h  and  0q  gives 

4 

3  +  5n  ’ 


l 


o  6  « 


8 


3  + 


80 


9  +  —  xy(n,a) 

vT 


(26) 


The  theoretical  relationship  between  fin  effectiveness 
€  and  the  non-dimensional  group  y  is  shown  in  figure  2.4 
for  the  cases  of  a  =  oo  and  1„ 
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FIG.  2.4  Fin  Effectiveness 


CHAPTER  III 


EXPERIMENTAL  WORK 


3.1  DESIGN  OF  APPARATUS 


In  view  of  the  approximations  used  in  the  theoretical 
section  of  this  thesis,  it  was  deemed  necessary  to  perform 
an  experimental  investigation  to  verify  them  for  the  con¬ 
ditions  considered.  From  the  analysis  of  Chapter  II  it 
was  predicted  that  similarity  solutions  to  the  governing 
boundary  layer  equations  are  possible  provided  a  power  law 
surface  temperature  exists.  The  corresponding  conduction 
problem  also  revealed  that  a  constant  heat  flux  at  the  fin 
root  would  provide  this  temperature  distribution.  There¬ 
fore  the  primary  purpose  of  the  test  apparatus  was  to  ob- 

I 

i 

tain  a  power  law  surface  temperature  on  a  vertical  fin 
immersed  in  an  isothermal,  quiescent  fluid. 

The  diagrammatic  and  pictorial  views  in  figures  3.1 
and  3.2  illustrate  the  test  apparatus  used.  It  consisted 
of  a  small  (5"  x  6-1/2"  x  8")  rectangular  aluminum  con¬ 
tainer  with  one  face  being  transparent  plastic  (PERSPEX) 

for  visual  observation.  A  tapered  aluminum  fin  was  in- 

* 

serted  in  a  vertical  position  through  the  top  of  the  container. 
The  three  fins  were  produced  from  4"  x  1"  aluminum 


plate  with  the  various  slenderness  ratios  obtained  by 
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FIG.  3.1  Diagrammatic  View  of  Apparatus 
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varying  the  fin  lengths.  A  straight  taper,  i.e.,  shape 
index  m  =  1,  was  used  for  each  fin  because  of  the  ease 
of  machining.  Also, since  the  theoretical  results  are 

I 

independent  of  shape  index  m,  any  shape  may  have  been 
chosen.  The  fins  varied  from  a  relatively  slender 
(w^Sr,  =  14,22)  to  a  rather  stubby  (w/ 6_  =  6,15)  shape 
to  reveal  the  effect  of  slenderness  ratio  on  the  problem. 

A  depth  of  one  inch  (figure  3,1)  at  the  root  of  each  fin 
was  not  tapered,  to  allow  mounting  space  and  thus  allowing 
only  the  tapered  portion  to  be  immersed  in  the  fluid. 

The  root  heating  of  each  fin  was  achieved  by  cem¬ 
enting  strands  of  #31  (0.0089"  diameter)  Greening  60-15 
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high  resistance  (8*25  ohras/ft.)  nickel-chromium  wire  on 
the  root  as  shown  in  figure  3.1.  The  cement  used  for  this 
purpose  was  BUDD  HT  (GA-100)  high  temperature  ceramic 
cement.  This  type  of  heater  along  with  the  one  inch 
untapered  portion  at  the  fin  root  was  considered  to  give 
a  constant  uniform  heat  flux  at  the  beginning  of  the  fin 
taper .  The  heaters  used  had  a  capacity  of  approximately 
200  watts  at  llOv  A»C»  A  variac  was  used  to  supply  the 
power  to  the  heaters  so  as  to  obtain  a  wide  range  of  heater 
input.  To  reduce  heat  losses  and  practically  eliminate 
heat  transfer  from  the  fin  to  the  box,  the  fin  root  was 
completely  enclosed  with  asbestos  as  illustrated  in 
figure  3.1. 

Since  the  problem  under  consideration  is  steady 
state  it  was  necessary  to  have  a  method  of  cooling  the 
ascending  boundary  layers  to  obtain  an  equilibrium  con¬ 
dition.  It  was  thought  that  thermoelectric  methods  would 
give  fine  control  and  consequently  twelve  thermoelectric 
cooling  modules  were  used.  The  Frigistor  modules  used 
are  rated  at  22.2  watts  for  a  maximum  current  of  thirty 
amperes  when  the  hot  face  is  held  at  27°C  and  the  cold 
face  is  thermally  insulated.  There  is  an  electrically 
insulating  film  on  both  sides  to  avoid  electrical  short 
circuits . 


These  modules  were  mounted  on  the  aluminum  container 


. 
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on  either  side  of  the  fin  root  as  shown  in  figure  3.1. 

Due  to  space  limitations  it  was  necessary  to  pyramid  the 
modules  with  two  layers  of  three  modules  on  each  side. 

These  twelve  modules  were  then  connected  in  series  and  the 
power  supplied  by  a  Hewlett  Packard  D.C.  power  supply 
rated  at  forty-five  amperes  and  eighteen  volts.  Since 
the  modules  require  an  adequate  heat  sink  in  order  to 
cool  properly,  a  1/4"  aluminum  plate  was  placed  over  each 
layer  of  modules.  The  modules  operate  on  relatively  small 
temperature  differences  between  "hot"  and  "cold"  sides 
and  since  the  "hot"  side  cannot  be  higher  than  100°C  a 
blower  was  used  to  cool  the  aluminum  plates. 

3.2  INSTRUMENTATION 

From  the  theoretical  analysis,  all  results  depend 
upon  the  temperature  index  n  and  therefore  it  was  nec¬ 
essary  to  obtain  a  tolerably  accurate  value  of  this  ex¬ 
perimentally.  Temperature  readings  were  therefore  required 
at  several  positions  along  the  fin  surface.  Embedding 
several  evenly-spaced  thermocouples  along  one  face  of 
the  fin  accomplished  this.  The  grooves  were  cut  laterally 
across  the  surface  such  that  the  individual  thermocouple 
wires  were  led  out  along  isotherms.  The  depth  was  just 
sufficient  to  allow  the  thermocouples  to  be  flush  with  the 
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surface  when  bonded  in  place  with  epoxy  glue*  Honeywell 
Type-J  iron-constantan  (0.010  inch  diameter)  butt  welded 
thermocouples  were  used  throughout  the  test  apparatus. 

Temperature  profiles  in  the  boundary  layer  adjacent 
to  the  fin  surface  were  obtained  with  a  probe  thermo¬ 
couple.  For  convenience  and  ease  of  locating  the  probe 
near  the  fin  surface,  the  wires  at  the  measuring  junction 
were  brought  to  a  blunt  point.  In  comparing  readings 
with  this  probe  and  an  "isothermal  junction"  probe,  there 
were  no  appreciable  differences  found;  thus  conduction 
errors  appear  to  have  been  tolerably  small. 

The  remaining  thermocouples  were  used  to  obtain  a 
bulk  fluid  temperature.  Since  it  was  not  known  if  an 
isothermal  fluid  could  be  attained,  several  thermocouples 
were  located  throughout  the  fluid.  These  were  all  located 
approximately  1/2 from  the  inside  surface  of  the  container. 

All  of  the  above-mentioned  thermocouples  were  con¬ 
nected  to  a  switching  unit  whose  single  output  was  measured 
with  a  potentiometer.  To  reduce  errors  in  thermocouple 
readings,  an  ice  bath  was  used  as  the  reference  point. 

3.3  EXPERIMENTAL  METHOD 


The  theoretical  prediction  that  all  results  are  de¬ 
pendent  upon  one  non-dimensional  group,  x>  suggested  that 
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the  tests  be  designed  to  verify  this  over  a  range  of  con¬ 
ditions  of  practical  interest 0 


3.3-1  TESTING  PROGRAM 


The  experimental  work  was  carried  out  in  three  phases. 
The  first  investigations  were  carried  out  to  test  the 
validity  of  the  n  versus  x  relationship.  The  non-dimen¬ 
sional  group  x  consists  of  three  separate,  non-dimensional 
groups;  namely  the  conductivity  ratio ,  the  slenderness 
ratio  and  a  variable  Prandtl-Grashof  number.  Two  values 
of  the  conductivity  ratio  were  obtained  by  using  both 
water  and  glycerine  as  the  bulk  fluids §  this  would  also 
give  an  indication  of  any  independent  effect  of  Prandtl 
number.  The  effect  of  variations  in  slenderness  ratio 
was  studied  by  using  the  three  different  fins  previously 
mentioned.  Variations  in  the  Prandtl-Grashof  number , 
other  than  by  using  two  different  fluids,  were  achieved 
by  changing  the  root  temperature  of  the  fin. 

Temperature  profiles  in  the  boundary  layer  adjacent 
to  the  fin  surface  were  measured  to  verify  the  profiles 
obtained  in  Chapter  II.  Four  profiles  were  measured,  two 
each  in  water  and  glycerine  for  one  fin  (w/6p  ~  14.22) . 
These  were  obtained  at  various  positions  along  the  fin 
surface  in  order  to  check  the  validity  of  the  similarity 


transformations . 
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The  last  phase  of  the  experimental  analysis  was 
obtaining  photographs  of  the  boundary  layer  for  two 
ranges  of  the  temperature  index  n.  To  obtain  these  a 
fin  was  mounted  in  a  transparent  plastic  (PERSPEX)  tank 
which  was  then  placed  in  a  Schlieren  apparatus  modified 
to  operate  as  a  shadowgraph .  The  purpose  of  these  tests 
was  to  examine  the  shape  of  the  boundary  layer  and  heat 
flux  curve  for  n  <  0.2  and  n  >1. 

3.3-2  PRELIMINARY  TESTS 


In  the  early  stages  of  the  experimental  work  several 
trial  runs  were  made  to  check  the  operation  of  the  apparatus. 
From  these  it  was  found  that  very  slight  variations  in 
surface  thermocouple  readings  made  it  extremely  difficult 
to  establish  a  reliable  value  of  the  temperature  index  n. 

Due  to  the  method  used  in  manufacturing  the  thermocouples 
it  was  deemed  necessary  to  calibrate  them  all.  These  were 
calibrated  at  both  the  freezing  and  boiling  points  of  water 
and  only  those  with  errors  of  less  than  0.1 °F  were  used. 

During  the  calibration  of  the  thermocouples  it  was 
noted  that  very  small  temperature  gradients  across  the 
switching  unit  would  give  false  readings  in  the  thermo¬ 
couples.  This  difficulty  was  eliminated  by  covering  the 
switching  unit  with  a  cardboard  box  to  maintain  isothermal 
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conditions  around  the  unit. 

Preliminary  tests  performed  in  water  showed  very 
erratic  readings  in  most  thermocouples .  This  was  attributed 
to  electrical  conduction  through  the  water  and  was  elimin¬ 
ated  by  coating  the  immersed  portion  of  all  thermocouples 
with  a  clear  acrylic  spray, 

3.3-3  TEST  PROCEDURE 


Considering  first  the  procedure  used  in  obtaining  the 
experimental  n  versus  x  relationship,  the  following  steps 
were  taken.  With  the  fin,  container  and  fluid  initially 
at  room  temperature  all  thermocouple  readings  were  re¬ 
corded.  If  any  were  in  appreciable  disagreement  with  the 
others  and  it  was  thought  to  be  due  to  malfunctioning  of 
those  thermocouples  for  some  particular  reason,  they  were 
disregarded  for  the  remainder  of  the  test.  Power  was  then 
supplied  to  the  modules  which  were  allowed  to  operate  for 
a  short  period  of  time  to  cool  the  fluid  slightly  before 
the  fin  heater  was  turned  on.  This  method  was  found  to 
reduce  the  time  required  to  obtain  steady  state  conditions. 

After  the  heater  had  been  switched  to  a  prescribed 
power  setting  and  allowed  to  operate  for  10  -  15  minutes 
the  power  to  the  modules  was  regulated  until  there  was 
no  appreciable  change  in  the  thermocouple  readings  over  a 
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10  -  15  minute  interval „  Once  steady- state  conditions 
had  been  attained  all  thermocouple  readings  were  recorded? 
the  power  to  the  heater  was  then  changed  and  the  procedure 
repeated 0 

In  obtaining  boundary  layer  temperature  profiles 
there  was  only  one  additional  step  required  to  those  given 
above „  After  steady  state  conditions  had  been  reached 
and  thermocouple  readings  taken,  the  probe  was  positioned 
as  accurately  as  possible  at  the  fin  surface  at  a  particu¬ 
lar  distance  from  the  fin  tip.  Probe  readings  were  then 
taken  for  equal  increments  of  distance  proceeding  from 
the  fin  surface  to  the  outer  region  of  the  thermal  bound¬ 
ary  layer 0 

The  visual  tests  consisted  of  setting  up  the  fin 
and  plastic  tank  in  the  Schlieren  apparatus  and  applying 
the  appropriate  power  to  obtain  the  desired  value  of 
temperature  index  n„  As  before,  the  value  of  n  was  de¬ 
termined  from  the  fin  surface  temperatures  and  the  bulk 
fluid  temperature..  Photographs  were  then  taken  of  the 
corresponding  boundary  layers  and  heat  flux  curves  as 


seen  on  the  screen „ 
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CHAPTER  IV 


DISCUSSION  OF  RESULTS 


4 „ 1  HEAT  TRANSFER  AND  SKIN  FRICTION 


The  eigenvalues  as  shown  in  figure  2.1  were  obtained 
using  seven  terms  in  each  of  the  two  infinite  series.  To 
note  the  effect  of  the  number  of  terms  in  the  series,  the 
eigenvalues  were  also  obtained  using  4,  5  and  6  terms. 
Maximum  variation  in  the  results  for  all  four  cases  was 
found  to  be  less  than  5  percent.  Also,  variations  in 
the  values  for  six  and  seven  terms  were  of  the  order  of 
1%,  indicating  that  seven  terms  are  sufficient  for  practi¬ 
cal  purposes  and  give  more  accurate  results  than  in  simi¬ 


lar  work 
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The  predominant  effects  of  n  on  the  heat  transfer 
and  skin  friction  results,  for  the  range  of  n  considered, 
is  clearly  shown  in  figure  2.1.  It  is  evident  that  in¬ 
creasing  n  causes  a  steepening  of  the  temperature  profile 
and  similarly  it  suggests  a  monotonically  decreasing 
tendency  for  the  velocity  profile.  Whether  these  trends 
would  continue  for  larger  values  of  n  cannot  be  established 
from  the  results  shown. 

The  independent  effect  of  Prandtl  number  on  any 


particular  fin  temperature  distribution  appears  quite 
slight  since  figure  2.1  reveals  only  a  small  change  in  heat 
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transfer  between  a  —  1  and  oo  over  the  range  of  n  considered 
This  tends  to  suggest  that  interpolation  in  the  range 
1  <  o  oo  ,  i»e„,  for  non-metallic  liquids,  should  be 
accurate.  Tolerably  accurate  extrapolation  into  the  gaseous 
region,  a  X  1,  is  also  reasonable  to  expect.  It  therefore 
appears  that  the  independent  effect  of  Prandtl  number  not 
contained  in  the  variable  Prandtl-Grashof  number  is  quite 
small  over  a  large  range  of  Prandtl  numbers.  Lefevre's 


results 


11 


for  an  isothermal  plate  suggest  that  y (n,a) 


11 


,  7(0, oo)  =  0.711 


is  more  dependent  on  n  than  a.  From 
and  7(0,0)  =  0.849.  Therefore,  results  presented  here 
appear  to  apply  for  all  but  the  smallest  Prandtl  numbers, 
which  correspond  to  the  liquid  metals. 

The  neglect  of  the  dissipation  terms  is  justified  since 

_  ^ 

the  Ostrach  number,  0  ,  is  of  the  order  of  10 

*  s 7 


4 . 2  TEMPERATURE  INDEX 


Perhaps  the  single  most  significant  result  to  arise 
from  the  theoretical  analysis  is  the  n  versus  ^  relation¬ 
ship  given  by  equation  (24) .  From  this  it  is  seen  that 
the  temperature  index  is  a  function  of  three  non-dimension- 
al  groups,  namely,  k^/k  ,  w/6R,  pg  ©R  w  //c  (1+a)  as  well 
as  the  heat  transfer  rate,  7,  at  the  fin  surface.  The 
conductivity  ratio  reflects  the  effect  of  conduction  through 
an  interface  between  two  different  conducting  media.  The 
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slenderness  ratio  reveals  perhaps  one  of  the  most  sur¬ 
prising  results  -  that  the  fin  temperature  distribution 

depends  on  overall  geometry  as  opposed  to  the  precise 

3  2 

form  of  the  fin  profile .  The  group  f3g  Gw  //c  ( 1-ha)  re- 

K. 

suits  from  the  Lefevre-type  transformation  and  may  be 
taken  as  a  generalized  "driving  potential" 0  This  group, 
unlike  the  Grashof  or  Rayleigh  numbers,  completely  speci¬ 
fies  conditions  over  the  entire  Prandtl  number  range.  It 
2 

becomes  a  Gr  as  a  — ►O  and  aGr  as  a  — *-oo  ,  both  of  which 
are  known  to  be  the  correct  forms  of  the  "driving  potential" 
for  these  extreme  values  of  a. 

The  theoretical  n  -  x  relation  in  figure  2.3  for  the 
range  0  <  n  <  3o0  illustrates  the  effective  independence 
of  Prandtl  number  in  the  range  1  <  a  <  go  .  This  independ¬ 
ence  implies  that  the  results  will  at  least  be  accurate 
for  all  non-metallic  fluids.  Table  I  is  a  summary  of 
test  results  from  which  values  of  x  were  calculated.  The 
corresponding  experimental  values  of  the  temperature  in¬ 
dex  n  were  obtained  from  the  slope  of  a  log-log  plot  of 
surface  temperature  versus  dimensionless  distance  along 
the  fin  surface.  The  root  temperature  was  calculated 
by  extrapolation  to  the  root  according  to  these  values 
of  temperature  and  temperature  index.  A  sample  calcula¬ 
tion  of  the  conversion  of  thermocouple  readings  may  be 
found  in  Appendix  C. 


££ 


1  \  a  '  &  :p  I  d' 


1  *nr  £/  j  f.  -tt  70  ji.c  i  o  ,exl 


- 


"  '3:-  !  I  :=  i  I  ■  r  „  7  9.  1*  ■ ! 3  .  c  i  9300 


■ 


39 


Slenderness 

Ratio 

(w/6r) 

14.22 

10.63 

6.15 

Conductivity 

Ratio 

(Vks) 

. 00264 

.00125 

.00261 

.00124 

.00265 

.00124 

Non-dimensional 

Temperature 

/Pg  eR  w3  \  _6 

(- 1-1 - )  *  10 

\c  (1+a)  / 

166  .  5 

5.61 

74.3 

1.86 

14.1 

1.17 

171.5 

1.06 

52.8 

( 

3.0 

18.2 

0.95 

159 

1.33 

96 

1.09 

25 

0.184 

156.5 

1.06 

161 

3.83 

35 

0.57 

162.5 

2.43 

23.7 

4.15 

92.2 

— 

165.5 

- 

— 

— 

i 

— • 

TABLE  I 


SUMMARY  OF  TEST  RESULTS 
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FIG,  4,1  Comparison  of  n  -  x  Relation 
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Figure  401  illustrates  the  comparison  of  theory  and 
experimental  results  for  the  range  0  <  n  <  1.5.  With 
the  exception  of  two  points  the  experimental  results  all 
lie  within  the  +  10%  error  band  as  denoted  by  the  broken 
curves „  Since  the  scatter  of  experimental  points  is 
fairly  random,  there  appears  to  be  very  little  effect 
due  to  the  changes  in  Prandtl  number  or  slenderness  ratio. 

One  basic  assumption  made  in  the  theoretical  analy¬ 
sis  was  that  of  quasi-one-dimensional  conduction.  This 
assumption  is  normally  valid  for  the  case  of  large 
slenderness  ratios.  It  was  therefore  somewhat  surprising 
to  find  from  the  experimental  results  of  figure  4.1, 
that  the  relatively  small  slenderness  ratio  of  6.15  gave 
no  greater  deviation  from  theory  than  did  the  higher 
ratios.  This  is  mainly  attributed  to  the  very  high  ratio 
of  fin  conductivity  to  fluid  conductivity  (380:1  in  water 
and  800:1  in  glycerine),  thus  giving  negligible  lateral 
temperature  gradients  in  the  fin.  What  would  occur  if 
the  conductivity  ratio  were  much  smaller  remains  to  be 
demonstrated . 

Throughout  the  analysis  the  effects  of  both  inclina¬ 
tion  and  curvature  were  neglected  on  the  basis  of  a  large 
slenderness  ratio.  The  effects  of  curvature  will  not 
be  discussed  here  since  the  experimental  results  were 
restricted  to  flat  surfaces.  However,  the  experimental 
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analysis  did  involve  inclination „  If  taken  into  account, 

the  main  contribution  due  to  inclination  will  occur  in 

the  x-momentum  equation  in  which  the  buoyancy  term  is 

replaced  by  (3g  cos  A(T  -  T  ) ,  where  the  fin  half  angle, 

A,  is  obtained  from  the  expression  6  /w  “  tan  A0  The 

JK 

additional  term,  cos  A,  would  occur  throughout  the 
analysis  as  (cos  A)  '  .  Considering  the  smallest  slender- 

ness  ratio  (w/6R  =  6015)  tested,  A  is  9°  14'  and  (cos  A) 
0.9967,  thus  introducing  an  error  of  only  0.325%.  There¬ 
fore  it  appears,  as  the  experimental  results  demonstrate, 
that  for  practical  cases  the  effects  of  inclination  are 
negligible.  It  is  probable  that  the  assumption  of  quasi- 
one-dimensional  conduction  will  be  violated  before  inclina¬ 
tion  has  any  significant  effect. 

4.3  VELOCITY  AND  TEMPERATURE  PROFILES 

The  theoretical  velocity  and  temperature  profiles 
in  the  boundary  layer  as  obtained  in  Chapter  II  are 
shown  in  figure  2.2.  These  profiles  reveal  the  familiar 
characteristics  of  a  large  Prandtl  number  system  in 
which  the  momentum  boundary  layer  extends  beyond  the 
thermal  boundary  layer.  Due  to  the  boundary  conditions 
for  this  system  the  temperature  profile  is  contained 
between  the  heated  fin  surface  and  the  velocity  maximum. 
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Utilizing  the  Lefevre-type  transformation  as  a — ►oo 
results  in  the  neglect  of  the  inertia  terms  but  not  the 
corresponding  advective  terms  in  the  energy  equation. 
Therefore  the  profiles  represent  the  solution,  within 
the  thermal  boundary  layer,  of  the  complete  energy 
equation  together  with  a  reduced  momentum  equation. 

The  velocity  profiles  were  not  determined  beyond 
the  thermal  boundary  layer.  However,  it  is  possible 
to  determine  the  velocity  profile  in  the  outer  region 
by  considering  the  complete  momentum  equation  (without 
the  buoyancy  term) .  The  inner  boundary  condition  for 
this  outer  region  would  be  taken  as  the  velocity  at  the 
edge  of  the  inner  layer.  Since  the  heat  transfer  and 
skin  friction  results  are  not  affected  by  the  resulting 
"tail"  of  the  velocity  profile,  it  was  not  considered 
in  the  present  analysis. 

A  comparison  of  theoretical  and  experimental  temp¬ 
erature  profiles  for  n  <  1  in  glycerine  is  presented  in 
figure  4.2.  Near  the  fin  surface  where  temperature 
gradients  are  highest  the  scatter  is  small  and  the  ex¬ 
perimental.  points  are  consistently  below  the  theoretical 
curve.  Most  of  this  effect  can  be  attributed  to  conduc¬ 
tion  losses  caused  by  the  blunt  point  of  the  probe 
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FIG.  4.2  Comparison  of  Temperature  Profile  (Glycerine) 
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FIG.  4.3  Comparison  of  Temperature  Profile  (Water) 
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thermocouple.  As  the  temperature  gradient  decreases, 
the  experimental  points  move  closer  to  the  theoretical 
curve  until  they  lie  very  nearly  along  it  in  the  central 
portion  of  the  thermal  boundary  layer.  However  in  the 
outer  region  the  scatter  becomes  much  more  noticeable. 

It  is  in  this  region  that  the  temperature  differences 
become  much  smaller  and  more  difficult  to  measure,  thus 
accounting  for  most  of  the  scatter. 

Using  the  same  fin  similar  results  for  n  >  1  in 
water  are  shown  in  figure  4.3.  In  these  tests  the  temp¬ 
erature  gradients  were  less  than  those  in  glycerine  and 
consequently  the  conduction  error  near  the  fin  surface 
should  be  less,  as  confirmed  in  the  figure.  The  scatter 
is  again  greatest  in  the  outer  region  of  the  thermal 
boundary  layer  where  temperature  differences  are  very 
small.  A  very  significant  result  from  this  figure  is 
the  close  agreement  between  experiment  and  theory  des¬ 
pite  the  fact  that  the  latter  is  for  a  =  oo  . 

As  previously  mentioned  the  temperature  profiles 
were  obtained  for  different  positions  from  the  leading 
edge  to  check  the  validity  of  the  similarity  transfor¬ 
mation.  From  figures  4.2  and  4.3  the  results  for  dif¬ 
ferent  heights  appear  to  be  quite  distinct  in  the  central 
and  outer  portions  of  the  boundary  layer.  This  variation 
of  results  with  changing  distance  from  the  leading  edge, 
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if  genuine,  implies  that  the  similarity  variable  £  is 
incorrect .  However  as  seen  from  the  figures,  this  effect 
does  not  appear  to  be  genuine  because  the  points  nearer 
the  leading  edge  form  the  upper  curve  in  glycerine  and 
the  lower  curve  in  water.  This  reversal  of  results  is 
not  likely  to  be  caused  by  changes  in  the  Prandtl  number. 
It  is  therefore  concluded  that  the  separation  is  mainly 
due  to  scatter  caused  by  errors  in  measuring  the  small 
temperature  differences. 

Since  the  scatter  of  experimental  results  is  no¬ 
where  large,  figures  4.2  and  4.3  together  reveal  the  over¬ 
all  effect  of  n  on  the  temperature  profiles.  As  predicted 
from  the  theory,  an  increase  in  n  from  about  0.4  in 
figure  4.2  to  about  1.35  in  figure  4.3  results  in  a 
steepening  of  the  profile. 

4.4  FIN  EFFECTIVENESS 

Utilizing  the  definition  of  fin  effectiveness  € 
given  in  Chapter  II  resulted  in  figure  2.4,  giving  the 
relationship  between  €  and  x-  BY  definition,  the  effect¬ 
iveness  may  be  greater  or  less  than  unity  depending  on 
whether  the  average  heat  flux  is  greater  or  less  than 
the  root  heat  flux.  From  equation  (20)  the  surface  heat 
flux,  h9Q,  varies  as  x  - —  ^  1  .  Thus  when  0_<  n  <  *jr  the 
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heat  flux  decreases  monotonically  from  the  leading  edge 
so  that  e  >  1.  For  the  case  of  n  =  ~  the  surface  heat 

D 

flux  becomes  uniform  and  €  =  10  Lastly,  when  n  >  ~  the 
heat  flux  increases  monotonically  from  the  leading  edge 
and  consequently  e  <  1. 

The  results  for  a  =  1  and  oo  are  effectively  the 
same  and  thus  imply  that  they  will  be  accurate  for  all 
non-metallic  liquids,  if  not  all  liquids.  It  is  revealed 
from  figure  2.4  that  a  change  in  effectiveness  is  more 
dependent  on  a  change  in  x  than  a  separate  change  in  a 
(with  x  fixed) .  This  indicates  that  the  effect  of  a  is 
largely  contained  in  x»  This  fact  is  further  demonstrated 
by  considering  the  results  of  the  following  tables  A  and  B. 


0 

n  =l 

7 

X 

00 

0  0  802 

0.353 

1 

0.773 

0.365 

0 

(0.920) 

(0.308) 

0 

X 

n 

€ 

00 

0.353 

0.20 

1.0 

1 

0.365 

0.206 

0.992 

0 

0.308 

0.179 

1.027 

TABLE  A  TABLE  B 

Although  results  are  not  available  for  Y(n,0),  other 
than  y (0,0)  =  0.849,  it  is  reasonable  to  expect  that  the 
shape  of  the  7  versus  n  curve  for  a  =  0  will  approximately 


•'  *'  :  '  *  ■  c v 

.  sbiupiX  XX*  Jon  11  tebivpxX  oXXXfilam-ooa 
n  .  pri  if  i  s  t  i  y  n  t  *:  ,nsrfo  B  r.o  zh^bncqob 


49 


follow  those  for  a  =  1  and  oo ,  at  least  for  small  n. 

It  is  on  this  basis  that  a  value  y(-|,0)  =  0.920  was 
chosen.  The  7  were  then  used  to  determine  the  corres¬ 
ponding  values  of  x  from  equation  (24) ,  as  shown  in 
table  A.  The  independent  effect  of  a  is  seen  to  be 
relatively  small  as  there  is  a  difference  of  only  12.5% 
between  the  values  of  x  f°r  0  =  0  and  o  =  00  . 

However 9  as  would  occur  in  most  practical  cases,  a 
value  of  x  would  be  known  and  it  would  be  desired  to 
evaluate  the  corresponding  fin  effectiveness  €.  Table 
B  considers  the  values  of  x  from  table  A,  for  which  the 
actual  values  of  n  and  €  are  0.2  and  1.0  respectively 
and  from  these  are  found  corresponding  values  of  n  and 
€  by  using  the  a  =  00  results.  The  values  of  e  obtained 
in  this  manner  demonstrate  that  the  independent  effect 
of  a  on  the  fin  effectiveness  is  very  small.  Thus  the 
results  shown  in  tables  A  and  B,  though  not  necessarily 
conclusive  for  all  n,  tend  to  reinforce  the  conclusion 
that  a  Prandtl  number  effect,  other  than  that  in  Xt  is 
sufficiently  small  to  allow  the  use  of  the  present  heat 
transfer  results  for  all  Prandtl  numbers. 

4.5  BOUNDARY  LAYER  SHAPES 

Photographs  obtained  in  the  Schlieren  apparatus  for 
the  cases  of  n  <  0.2  and  n  >  1  are  shown  in  figures  4.4 
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and  4.5  respectively.  The  purpose  of  these  was  to  verify 
the  shape  of  the  boundary  layer  and  heat  flux  curve  for 
differing  values  of  n„ 

As  mentioned  in  the  previous  section,  the  heat  flux 
decreases  monotonically  from  the  leading  edge  when 
0  <  n  <  “|  and  increases  monotonically  when  n  >  These 
two  facts  are  observed  in  figures  404  and  4„5  respectively, 

An  interesting  feature  of  boundary  layer  growth  is 
revealed  by  examination  of  the  similarity  variable  £. 


n 


The  boundary  layer  thickness  is  found  to  vary  as  x  — 
so  that  when  n  <  1  the  thickness  increases  monotonically 
with  height  and  decreases  monotonically  when  n  >  1.  Of 
course,  with  n  =  0  the  familiar  isothermal  case  occurs  in 


which  the  thickness  varies  as  x 


1/4 


When  n  =  1  the 


thickness  is  independent  of  height  and  therefore  the 
boundary  layer  must  be  parallel  to  the  surface.  Visual 
confirmation  of  the  effect  of  n  on  the  boundary  layer 
growth  is  obtained  from  the  Schlieren  photographs  in 
figures  4,4  and  4.5. 


v  lib 


- 


FIG.  4.4 


Schlieren  Photograph  (n  <  0.2) 
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Schlieren  Photograph  (n  >  1) 


FIG.  4.5 


CHAPTER  V 


CONCLUSIONS  AND  RECOMMENDAT IONS 


5.1  CONCLUSIONS 


This  thesis  has  described  a  theoretical  and  experi¬ 
mental  investigation  of  the  combined  conduction-convection 
problem  of  heat  transfer  through  and  from  a  tapered,  down¬ 
ward-projecting  fin  immersed  in  a  quiescent,  isothermal 
fluid.  The  theoretical  analysis  has  demonstrated  that 
a  combined  problem  of  this  type  may  be  solved  "exactly" 
if  certain  simplifications  are  made  in  the  individual 
conduction  and  convection  problems.  A  power  law  solution 
for  surface  temperature  was  revealed  after  solving  the 
convection  problem  separately  and  then  applying  the  con¬ 
dition  of  continuity  of  temperature  and  heat  flux  at  the 
common  interface. 

The  overall  heat  transfer  problem  was  found  to  be 
governed  by  a  single  non-dimensional  group  which  con¬ 
sists  of  a  slenderness  ratio,  conductivity  ratio  and  a 
variable  Prandtl-Grashof  number.  The  results  were  thus 
found  to  be  dependent  only  on  the  overall  geometry  as 
opposed  to  the  exact  fin  profile.  The  separate  effect 
of  Prandtl  number  on  the  results  was  found  to  be  vir¬ 
tually  negligible.  In  fact  the  temperature  index  and 
fin  effectiveness  were  much  more  sensitive  to  changes  in 


v 
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root  temperature ,  fin  material  and  overall  geometry 
than  a  separate  change  in  the  surrounding  fluid 0  This 
separate  independence  of  Prandtl  number  applies  over  the 
range  1  <  a  <L  oo  ,  if  not  the  entire  range „ 

The  majority  of  the  theoretical  predictions  were 
well  borne  out  by  the  experimental  results  obtained.  A 
comparison  of  theoretical  and  experimental  results  for 
the  n  versus  x  relationship  revealed  agreement  to  within 
10%  or  better o  The  experimental  results  clearly  verified 
the  overall  trend. 

Comparison  of  experimental  and  theoretical  tempera¬ 
ture  profiles  for  water  and  glycerine  revealed  good 
overall  agreement  despite  the  fact  the  theoretical  pro¬ 
files  were  for  o  =  oo.  Discrepancies,  although  never 
large,  were  greatest  in  the  outer  region  of  the  thermal 
boundary  layer  where  temperature  differences  were  very 
small  and  thus  difficult  to  measure.  Since  the  tempera¬ 
ture  profiles  were  found  to  be  independent  of  distance 
from  the  leading  edge  of  the  fin,  the  validity  of  the 
similarity  variable  £  was  revealed. 

The  results  and  hence  the  conclusions,  are  only  valid 
for  the  case  of  a  single  vertical  fin  or  banks  of  vertical 
fins  in  which  there  are  negligible  effects  of  one  fin  upon 
another.  The  results  for  relatively  small  slenderness 
ratios  are  only  valid  where  the  fin  conductivity  is 
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sufficiently  greater  than  that  of  the  surrounding  fluid, 
as  in  the  experiments  reported  here  and  in  most  practical 
applications . 

5.2  RECOMMENDATIONS 


The  analysis  and  thus  the  results,  are  limited  to 
the  behavior  of  a  single  vertical  fin.  A  logical  exten¬ 
sion  would  be  to  consider  multi-fin  systems  which  are 
more  likely  to  occur  in  practice.  In  this  case  an 
additional  parameter  enters  into  the  problem,  namely, 
the  pitch/height  ratio?  this  gives  an  indication  of  the 
effect  of  mutual  proximity. 

The  only  orientation  considered  for  the  fin  in  this 
analysis  was  that  in  a  vertical  downward-projecting  posi¬ 
tion.  Other  positions  such  as  horizontal,  upward-project¬ 
ing  or  angles  in  between  these,  could  be  studied.  Also 
since  no  results  were  obtained  for  fins  having  curvature, 
studies  could  be  carried  out  to  determine  these  effects 
which  would  probably  be  some  function  of  the  radius  of 
curvature.  A  cooled  fin  for  which  the  flow  would  be 
reversed,  at  least  for  the  vertical  position,  is  another 
possibility  for  further  analysis. 

The  validity  of  assumptions  made  in  the  theoretical 
section  was  ascertained  with  the  use  of  experimental 
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temperature  indexes  and  profiles  only*  A  study  of  the 
experimental  velocity  profile  would  further  verify  the 
results  and  may  also  indicate  the  effect  of  neglecting 
the  inertia  terms  for  liquids  whose  Prandtl  number  is 
small.  However ,  as  in  most  free  convection  problems  the 
velocities  are  small  and  very  difficult  to  measure  direct¬ 
ly.  One  feasible  method  would  likely  be  the  use  of  the 
fibre  anemometer „  Dye  streaks  in  glycerine  and  suspended 
particles  in  water  will  give  indications  of  bulk  fluid 
velocities . 

The  experimental  results  revealing  the  effects  of 
slenderness  ratio  on  the  problem  were  all  obtained  for 
the  case  of  a  very  high  ratio  of  fin  to  fluid  conductivity. 
Further  tests  with  much  smaller  ratios  should  be  per¬ 
formed  to  determine  the  possible  effects  of  lateral  temp¬ 
erature  gradients  in  the  fin. 

Further  investigations  could  be  performed  to  deter¬ 
mine  whether  or  not  the  effect  of  Prandtl  number  is  suf¬ 
ficiently  contained  in  the  variable  Prandtl-Grashof  number. 
A  series  of  tests  could  therefore  be  performed  in  air  and 
liquid  metals  to  give  an  indication  as  to  whether  or  not 
the  predictions  discussed  in  chapter  IV  are  valid* 
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APPENDIX  A 


FORTRAN  IV  PROGRAM  FOR  EIGENVALUE  SOLUTION 

The  computer  program  used  to  compute  the  values  of 
a  and  7  versus  n  follows  this  introduction „  Assuming 
initial  numerical  values  for  a(=A)  and  7(=B)  the  pro¬ 
gram  sets  n(=-H)  to  zero  and  calculates  the  values  of 


the  coefficients  ar  (=Ar)  ,  b  (=Br  1),  b^C^Br  2), 

c  (  =  Cr)  ,  d  ( =  Ds  1)  and  g  Ds  2)  ,  From  the  d  and 
r  s  s  s 


are  formed  the  products 


which  are  then  transferred  to  the  subroutine  (EULER)  to 
check  the  convergence  and  apply  Euler’s  transformation  if 
necessary„  This  same  procedure  is  carried  out  for  the 
products 


g 


s 


The  series 


are  both  summed  and  checked  against  their  known  and  speci¬ 
fied  correct  value „  If  they  are  incorrect  either  a  or  7 
is  incremented  and  the  coefficients  and  series  are  then 
recalculated,,  This  iteration  procedure  is  repeated  until 
both  S ( I)  and  R(J)  assume  acceptable  values 0  Once  the 
latter  are  obtained  n(=H)  is  incremented  a  given  amount 
and  the  whole  procedure  is  repeated s  using  as  initial 
values  of  a  and  7  the  values  obtained  for  the  previous  n0 
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1C0)**3 ) /81 •) )/( 3**C0* *1.333 3 3733 ) 
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T-  3  )  /9  •  )/7*+BOl*(-6.*C4/(R.*C0)+70.*(  (C7/CO)**?  )/9.-770.*(  (ri**7)*C7 
2/(CC**3)  )/27.+770.*( (Cl/C0)**4)/243.)  )/(3.*C0**1.66666667) 

D51  =  (B51 /120.-B41*Cl/( ] 2 . *CO ) +B3 1  * ( -2 . *C2 /CO  +  R . * (  (Cl/CO)**7)  ) /6 • +R 
l71*(6**Cl*C7/( CO** 7 )— 4«*(Cl/C0)**R)/7*+R0l*(— 2**C3/C0  +  6.*C1*C4/(C0 
7**2 )~12  **C1*( C2**2 ) / ( CO**3 ) +20**C2* ( Cl**3 ) / ( CO**4 ) -6** ( Cl /CO )**5 )  ) 
3 / ( 3 • *CO**2 ) 

O61=(R61/720.“7.*B3]*Cl/(^60.*C0)+P41*(-7**C2/(8.*C0)+33.*((Cl/CO) 
1 **  7  )  /9 •  )/74.+B71*(70.*Cl*C7/(Q**C0**7)-466.*(  (Cl/CO)**3)/pl.  ) /6.+R 
?71*(-7**C4/ ( ** • *C0 ) +RP .* (  (C7/C0)**7) /Q.-466.*C7*(Cl**7) /(27**CO**3) 
3  +  1  820.  *  (  (Cl  /CO  )  **4  )  /74^».  )  /  2  .  +  R0  1  *  ( -7  .  *f6  /  (  7  •  *C0  )  +70  •  *  (  Cl  [*C5+C?  *C4  ) 
4/(9.*C0**7)-455.*(‘2.*C4*(Cl**7)  +  (C7**7  )  )/(81  • *C0**^  )+^460.*(Cl**7) 
^* ( C?**7 )  /  (pi  .*r0**4)-'*4*P0**C7*(C1  **4)/(779.*C0**3)+76076.*(  (Cl  /  CO 
6 ) **6 ) /6561 • ) )/(3.*C0* *2*33333333) 

D02=B02/ (3.*C0**.33333337 ) 

D12=(B12-2.*B02*C1 /(3.*C0) ) / ( 7 . *C0**. 66666666 ) 

D27=(B02*( (Cl /C0)**7-C7/C0)-B1 7*C1 /CO) /(7.*C0) 

0?7  =  (B37/6*+Rl  7*  (~4**C7  /  (  7**C0  )  +1  4**(  (  Cl  /CO  )**?  )  /Q.  )  4-R07*  (  78  **C1  *C 
17/(9* *C0**2 ) “1 40 .* ( (Cl/C0)**3)/81. ) ) / ( 3 . *C0** 1 . 3 3 3 7 3 3 3 3 ) 

D4?= ( R4  7 /74*-3  **B37*C1 / ( 1 8  **C0 )+R17*(^0.*Cl*C7/(o.*C0**7)-770.*((C 
71  /CO)**^)/pi  #  )+B07*(~c’**C4/(7.*C0)+70.*(  (C7/CO)**7)/q#-22o.*C7*(Cl 
3**2 ) / ( 27  **C0**3 ) +7  70  *  * (  (Cl/C0)**4)/748.)  )/(3**C0**l  .66666667) 
D67=(B32/170*~B43*C1 /(17**C0)+B37*(-7**C7/C0+P.*( (Cl/C0)**7) )/6*+B 
1 1 7*( ~7**C4/C0+8**( (C7/C0)**2)-12.*C7*(Cl**7)/(CO**7)+5.*( ( Cl /CO ) ** 
74 )  ) +B07* ( -7  **C5 /CO+6*  *C1 *C4/ ( CO** 7 ) -1 7  **C1 * (C7**7 ) / ( C0**3 ) +70. *C7* 
3  (  0**3  )  /  (C0**4  )  -6**  (  (  Cl  /  CO  )  **S  )  )  )/(3.*C0**7) 

067=(B67/770.“7.*B67*C1/(760.*C0)+B47*(-7.*C7/(8.*C0)+3^.*((C1/C0) 
l**7)/9. ) /74.+B32* ( 70. *C1*C2/ ( 9.*C0**7 )-435**( ( Cl /CO ) **3 ) /fi 1 . )/6.+R 
217*(-7.*C5/(7.*C0)+70.*C1*C4/(Q.*C0**7)-4c;6**C1*  (  C7**7  )  /  (  77.  *C  0**7 
3 )+7?80.*C7* (Cl **3 ) / ( ? 47 .*C0**4 )~691 6**( (Cl /CO ) **6 ) /729. )+B07* (-7.* 

4C6/(3**C0)+70**(Cl*C6+C2*C4)/(9.*CO**2)-46S.*(7.*C4*(Cl**7)+(C3**7 
c'))/(81**C0**7)+54^0.*(Cl**7)*(C2**7)/(81.*C0**4)--746P0.*C7*(Cl**4) 
6/ ( 729  «  *C0**5 ) +76076  *  * ( (Cl/C0)**6)/6361* )  )  / ( 3 • *C 0**2*73333333) 

RO 1 =D0 1 *G0 
R11=D 11*61 
R21=D21*G2 
R31 =D31 *63 
R4 1 =D4 1*64 
R5 1=D5 1 *65 
R6 1 =D6 1*66 
T0=R01 
T  1  =R  1 1 
T  2  =  R21 
T  7  =R  3 1 
T4=R41 
T5=R51 
T6=R61 

CALL  EULER(T0*T1»T2*T3*T4*T5*T6) 

R01 =T0 
R 1 1 =T 1 
R2 1 =T2 
R31 =T3 
R41 =T4 
R5 1 =T5 
R6 1 =T6 
R07=D07*60 
R12=D17*61 
R??=D22*G2 
R72=D32*63 
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R42=D42*64 

R52=D52*G5 

R62=D62*G6 
TO  =R02 
T 1 =R 1 2 
T2=R22 
T3=P32 
T4=R42 
T5  =R52 
T6=R62 

CALL  EULER ( TOtTl *T2 ,T3*T4*T5 *T6 ) 

R02=T0 

R 1 2=T 1 

R22=T2 

R32=T3 

R42=T4 

R52=T5 

R62=T6 

S( I )=R01+R11+R21+R31+R41+R51+R61 
R  (  J  )  =R02+R12+R22+R32+R42+R52+R62 
IF (D.GT.0.0 )GO  TO  40 
IF(M.GT.O)  GO  TO  35 
IF(N.GT.O)  GO  TO  30 
U  =  0.2*< ABS(S< I )  —  1 • 0  )  ) 

A=  A+G*U*  A 
N  =  N  4- 1 
1=1  +  1 
J  =  J  +  1 
GO  TO  20 

30  Z  =  APS  (  S  (  I  )  -1 • 0 ) 

W*  ABS ( S ( I -1 )  —  1 • 0  > 

IF (Z.GT.W)G=-G 

IF (Z.LE.O.ODl )  GO  TO  46 
45  U=0.2*< ABS<S< I >-1.0> > 

A=A+2.*G*U*A 
M  =  M+ 1 
J  =  J+1 
1  =  1  +  1 
N  =  0 

GO  TO  20 
35  C=C+1*0 

IF(L*GT.O)GO  TO  31 
V=0.2*( ABS(  ( R ( J ) +A ) / A )  ) 
r=r+T*V*B 
1=1  +  1 
J=  J+l 
L  =  L  + 1 
GO  TO  20 

31  0=  ABS ( R ( J ) + A ) 

X=ABS(R( J-l >+A> 

I F ( 0*GT • X ) T  =— T 

IF (Q.LE.0.001 )  GO  TO  72 
73  V=0.2*( ARS(  (R( J)+A) /A)  ) 

R  =  B+2  •  *T *V*B 
1=1+1 
J  =  J+1 
L  =  0 

D=D+1.0 
GO  TO  20 
M  =  0 
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IF ( RR.GT.0.0 )  GO  TO  41 
IF ( SS.GT.0.0 )  GO  TO  42 
U=  0. 2* ( ABS  ( (R( J )+A> /A  >  ) 

A= A+G*U* A 
SS*SS+1 .0 
1  =  1  +  1 
J= J+l 
GO  TO  20 
0=  ARS ( R ( J ) +A ) 

X= ABS ( R ( J-l ) + A ) 

IF (O.GT .X )G=-G 
I F ( Q. LF • 0*  001 )  GO  TO  70 
U=0.2*( ABS ( ( R ( J ) +A ) / A ) ) 

A=  A  +  2  .*G*U* A 
RR=RR+1 .0 
1=1+1 
J  =  J+1 
GO  TO  20 

ss=o.o 

IF(F.GT.0*0)  GO  TO  43 
V=  0»  2* ( ABS ( S (  I )-l .0)  ) 

B=p+T*V*R 
1  =  1  +  1 
J  =  J+1 
E=F+1 *0 
GO  TO  20 
Z= ABS ( S ( I )-1.0) 

W= ABS ( S ( 1-1 )-l .0 ) 

IF(Z.GT.W)T=-T 

IF ( Z*LE* 0.001 )  GO  TO  74 

V=0.2*< ABS(S( I ) — 1 • 0 ) ) 

B=B+2.*T  *V*B 
1  =  1  +  1 
J  =  J+1 
F=0. 0 
RR  =0 • 0 
D=  0  •  0 
GO  TO  20 
0=  ABS ( R ( J ) +A ) 

IFfO.LE. 0*001  )  GO  TO  10 
GO  TO  45 
Z  =  ABS ( S (  I )-1.0) 

IF ( Z.LE.0.001 )  GO  TO  10 
GO  TO  73 
Z  =  ABS ( S ( I ) —1  *  0 ) 

IF (Z.LF. 0.001 )  GO  TO  10 
GO  TO  71 
Q=  ABS ( R ( J ) +A ) 

IF (O.LF. 0.001  )  GO  TO  10 
GO  TO  75 

WRTTE(6*52)R01*R11*R?1,R31*R41*R51*R61 
WRITE(6*52)R02*R12*R22*R32*R42*R52>R62 
WRITE  (6*51 )A*R*S( I > *R( J) 

CONTINUE 

STOP 

FORMAT ( 2X .23HTEMPFRATURE  EXPONENT  =  *F8.4/) 

FORMAT ( 2  X  *  9H  ALPHA  =  *F12.8*9H  GAMMA  =  *F12.8*8H  SUM1  =  F12.8*8H  S 
1 UM  2  =  *  F  1  2  •  8  /  ) 

FORMAT (5X,6F12.6*F14.8) 

END 
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SIRFTC  EULER  NOL I  FT  *  NODECK 

FUBROIJT  INF  EULER <  TO *T1*T?*T3  *T4*T5  *T6 ) 

DIMENSION  RO ( 250 ) *Rl (?F0)>R2(250)9R*(?60),R4(250)*R5(250> *R6 ( 7^0 ) 
AB  =  1  . 

MN  =  1 


GO 

TO 

500 

GO 

TO 

500 

GO 

TO 

500 

GO 

TO 

500 

GO 

TO 

500 

GO 

TO 

500 

R0(mn)=T0 
R 1  ( MN ) =T 1 
R2 ( MN )  =  T2 
RF (MN )  =T3 
R4 (MN ) =T4 
R  5 (MN )=T5 
R6 (MN  )=T6 

IF(  ABS(RKMN)  )  .  GF  .  ABS  (  RO  (  MN  )  )  ) 

IF(ABS(R2(MN) ) .GF.ABS(R1 (MN) ) ) 

IF ( ABS ( R3(MN ) ) • GF . ABS ( R2 ( MN ) ) ) 

I F ( ABS ( R4 ( MN ) ) • GF . ABS ( R3 ( MN ) ) ) 

IF  (  ABS ( R5 (MN  >  ) .GF . ABS ( R4 ( MN  >  )  ) 

I F ( ABS ( R6 ( MN ) ) .GE.ABS(R5(MN) > ) 

NM  =  MN 
GO  TO  501 
5  00  MN  =MN+ 1 
NM  =  50 
JK  =  50 

R6 ( NM ) = • 500 
600  F0=R0 (MN-1 ) 

F] =R0(MN-1 )+Rl (MN-1 ) 

F2=F1+R 1 (MN-1 ) +R2 (MN-1 > 

F3=F2+R1 (MN-1 ) +2.*R2( MN-1 )+R3( MN-1 ) 

F4=F3+R 1 (MN-1 )+3.*R2(MN-l )+3«*R3(MN-l )+R4(MN-l ) 

F5=F4+R1 (MN-1 )+4.*R2(MN-l )+6.*R3 (MN-1 )+4«*R4(MN-l )+R5(MN-l ) 

F6  =  F  5+R 1  (MN-1 ) +5 .*R? ( MN-1  ) +10.*R? ( MN-1 ) +10.*R4(MN-1 ) +5 • *R R { MN-1 ) +R 
16 ( MN-1 ) 

RO ( MN ) =F0/2 • 

R1  (MN )=Fl/4. 

R2 (MN)=F2/8. 

R3 (MN )=F3/16. 

R4 (MN ) =F4/32 • 

R5 (MN )=F5/64. 

R6 (MN )=F6/128. 

I F ( ABS ( R 1 ( MN ) ) • GT • ABS ( RO ( MN ) ) ) 

I F ( ABS ( R  2 ( MN ) ) .GT.AB$(R1 (MN)  ) ) 

I F ( ABS ( R  3 ( MN )  ) .GE • ABS ( R2 ( MN )  )  ) 

I F ( ABS ( R4 ( MN ) > .GF. ABS ( R3 ( MN ) ) ) 

I F ( ABS ( R5 ( MN ) ) .GF. APS <P4 ( MN > ) > 

IF(ABS(R6(MN) ) .GE. ABS ( R 5 ( MN > )) 

I F ( ABS ( R6 ( MN ) ) . LT . ABS ( R6 ( NM ) ) )NM=MN 
401  CONTINUE 
MN=MN+1 


GO 

TO 

401 

GO 

TO 

401 

GO 

TO 

401 

GO 

TO 

401 

GO 

TO 

401 

GO 
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R6 l MN ) =F6/64 I 
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)  NM=l 

MN 

CONTINUE 
MN=MN+1 
KK=MN-3 
F2=R2(KK  ) 

F^>=R2  (KK  )+R3  (  KK  > 

F4  =  F3+R3 (KK )  +  R  4  (  K  K  ) 

F5  =  F4+R  3 (KK)+2.*R4(KK>+R5(KK) 

F6  =  F  5+R  3 (KK )+3 .*R4 ( KK  > +3 • *R 5 ( KK ) +R6 (KK  ) 
RO ( MN ) =R0 ( KK ) 

R 1 ( MN ) =R 1 ( KK ) 

R2 (MN )=F2/2. 

R3 (MN)=F3/4. 

R4 ( MN ) =F4/8 • 

R5 (MN ) =F5/16. 

R6 ( MN ) =F6  /32 • 

IF ( ABS ( R1 ( MN ) ) .GE. ABS ( RO ( MN )  )  ) 
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IF ( ABS ( R3 (MN ) ) • GE • ABS ( R2 ( MN ) ) ) 
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IF ( ABS ( R5 (MN > > .GE. ABS ( R4 ( MN > ) ) 

I F ( ABS ( R6 ( MN )  ) .GF • ABS ( R5 ( MN )  )  ) 

IF ( ARS( R6 (MN ) ) • L T . ARS ( R6 ( NM ) ) ) NM=MN 
CONT INUF 
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F3=R3 (KK  ) 

F4  =  R3 (KK )+R4( KK  ) 

FS=F4+R4 (KK )+R5 ( KK  ) 

F6  =  F?  +  R4(KK )+2.*R5 ( KK ) +R6 ( KK  ) 

RO(MN ) =R0 ( KK ) 

R 1 (MN)=R1 (KK) 

R2 ( MN ) =R2 ( KK ) 

R3 (MN)=F3/2. 

R4 (MN ) =F4/4. 

R5 (MN)=F5/B. 

R6 ( MN ) =  F6 / 1 6 • 

IF  (  ARS ( Rl (MN )  ) .GF. APS (RO(Mf 
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R4 (MN ) =F4/2 . 

R5 (MN ) =F5/4. 

R6 (MN)=F6/P. 

IF ( ARS( R1 (MN )  ) .GF.APS (R0(MN  > 
IF ( ABS ( R2(MN ) ) .GF. ABS ( R1 (MN ) 
I F ( ABS (  R3  ( MN )  ) .GF. ABS ( R2 ( MN ) 
IF(ABS(R4(MN) ) .GE.ABS(R3(MN) 
IF ( ABS ( R5 ( MN ) ) . GF • ABS ( R4 ( MN > 
IE(ARS(R6(MN)).GE.ARS(R5(MN) 
IF ( ABS( R6(MN> > • LT . ABS ( R6 ( NM > 
CONTINUE 
MN=MN+1 
F5=R5(KK ) 

F6  =  R5 IKK  )+R6( KK > 

RO ( MN ) =R0 ( KK ) 

R  1  ( MN ) =R 1 ( KK ) 

R2 ( MN ) =P2 ( KK ) 

R3 (MN ) =  R3  ( KK ) 

R4 (MN)=R4(KK) 

R5 (mn)=F5/2. 

R6 (mn ) =F6/4. 

IF (ABS( R1 (MN > > .GE.ABS (R0(MN ) 
IF ( ABS ( R2 (MN ) ) .GE.ABS (R1 (MN) 
I F ( ABS ( R3 ( MN ) ) .GE . ABS ( R2 ( MN ) 
I F ( ABS ( R4 ( MN )  ) .GF . ABS ( R3 ( MN ) 
I  F ( ABS ( R  5 ( MN )  ) . GF • ABS ( R4 ( MN ) 
I F ( ABS ( R6 ( MN ) ) • GF • ABS ( R 5 ( MN ) 
IF ( ABS ( R6 ( MN )  )  .  LT . ABS ( R6 ( NM ) 
CONTINIJF 

IF(NM.NF.JK)  GO  TO  501 
MN  =KK  +  1 
GO  TO  500 
TO=RO ( NM ) 

T 1=R 1 ( NM ) 

T2=R2 ( NM ) 

T3=R3 (NM ) 

T4=R4( NM ) 

T5=R5 ( NM ) 

T  6  =  R  6  (  N  M  ) 

RETURN 

END 
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APPENDIX  B 

FORTRAN  IV  PROGRAM  FOR  VELOCITY  AND  TEMPERATURE  PROFILES 


Once  the  eigenvalues  a  and  y  have  been  obtained  it 
is  a  simple  matter  to  calculate  the  velocity  and  tempera- 
ture  profiles  by  forward  integration „  The  computer  program 
used  for  this  purpose  follows  this  brief  introduction. 

The  values  of  the  eigenvalues  a(=C(K))  and  y(  =  D(K)) 
are  first  read  in  for  each  value  of  n(  =  K)„  Then  utilizing 
a  fourth  degree  polynomial  the  program  calculates  corres¬ 
ponding  values  for  the  following,  versus  £(==XX): 

F(i)  =  FO  ( I) 

F  0  (^)  =.  Fl  ( I) 

F  0  '  (?)  =  F2  (I) 

F8  8  '  H)  =  F3  ( I) 

FI¥(£)  =  F4(I) 


V 


. 
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SI8FTC  VFLTEM  NOL I  ST  * NODEC K 

DIMENSION  C(300)*D(20n),FO(300)*Fl(7nn)*F2(30n),F^(?nn),F4(300)*Fs 
1 (SCO) 

N  =  31 
E=0.^5 

READ (5 *51 > ( C ( K  > *D(K> *K  =  1 *N> 

DO  10  K  =  1  *N 
P  =  K 

H* (P-i . ) /in. 

WR I TE ( 6  *  50 ) H 
WRITE (6 *52) 

A  =  C  (  K ) 

R=D ( K ) 

20  A  2  = A 

A  5  =  —  1 .0 
A4  =  R 

A6  =  -4 •  *H*A 

A7=4 • *H  + (  7 • *H-3 • )*B*A 
A8  = (3.-15. *H)*B 

A9=( 1 8.-10. *H)*(4.*H*A*A>+( 15.*H-3. >*R*B 

A10  =  176**H*H*A-240.*H*A+70.*H*H*A*A*R-‘7  40.*H*A*A*P  +  Q0.*A*A*R 
A  1  !=1446.*H*A*R-795.*H*H*A*R-715.*A*R-176.*H*H+240.*H 
A  1 2=-2R0 •  * ( H**3 )  *(  A **2 ) +2024 .*H*H* ( A**7  ) +1 0 2 Q . *H*H* A*R*R- 2 1 1 8.*H*A 
1*R*R-216  6.*H*R+12;>3.*h*H*A-45  36.*H* ( A **7 ) +5 i 3 . * a*E*R+2 1 5 • *B 
DO  5  1=1 *4 
TT  =  I 

X  I  =  ( TT— 1 • ) /20 • 

F0(I)=A2*(XT**2)/2.+A2*(XI**2) /6  • +  A  4* ( X  I **4 ) / 24 . +A 6* ( X  I **6 ) /720.  +  A 
17* (X 1**7) /5040.+AR* (XI**R ) / 40320. 

Fl  ( I ) =A2*X I+A3*( X 1**2 ) /2.+A4* ( XI** 2 ) /6 .+A6* ( X 1**5 ) /120.+A7* ( X 1**6 ) 
l/720.+AR*(XI**7) /5040.+A9* ( X 1**8 ) /40220. 

F2 ( I )=A2+A3*X I +A4* ( X  I **2 ) /2 .+A6* ( X 1**4 ) /24 • +  A7* ( X  I **5 ) / 1 2  0 .+A8* ( X  I 
1  **6  ) /7  20.+A9* ( X 1**7  )  /8040.+A 1 0*  (  XI**R  )  /403?0. 

F3 ( I ) =A3+A4*X I+A6* ( X  1**2 ) /6 .+A7* ( X 1**4 )/24.+AR*(XI***)/120.+A9*(XI 
1 **6 ) /720.+A1 0* ( X 1**7 ) /5040.+A1 1* ( X I**R ) /4P27H . 

F4( I )=A4+A6*(XI**7 ) /? • +A7* ( X T **3 ) /6 . +A 8* ( X  I **4 ) / 24 .  +  A9* ( X  I ** 5  )  / 1 20 
1 *+Al 0* ( X 1**6 )/720.+All*(XI**7)/5040.+A12*(XI**8)/40?20* 

F5 ( I ) =4 • *H*F3 ( I  ) *F 1 (  I  )-(H+2. )*F4(  I )*F0(  I  > 

F  WRITF  (6*F7 )XI *F0( I ) *F1  ( I )  *F?( I  )  ,F2 ( I )  ,F4( I )  *FF(  I ) 

1=5 

15  11=1-1 

12=1-2 
13=1-3 
14=1-4 
TT  =  I 

X I  =  ( TT— 1 • ) /20« 

F0(  I )=F0( II  )+F*(Fl  (  n  )*0.5*(F1  (  II  )-Fl  (  12)  )  +  5.*(Fl  (  II  )-2.*Fl  (  I2)+Fl 
1 (  I3))/12,+3.*(F1 ( II )-3.*Fl ( I2)+3.*F1 ( 13 )“F1 ( 14) ) /8. ) 

Fl ( I)=F1 (II )+E*(F2( II )+0.5*(F2( II )“F2( 12) )+5.*(F2( II )“2.*F2( I2)+F2 
l(l3))/12.+3.*(F2(Il)-3.*F2(l2)+3.*F2(I3)-F2(l4))/8*) 

F2 (I )=F2 ( H)+E*(F3 ( II )+0«5*( FR(  1 1 )-F3( 12 ) )+5.*(F3( II  )-2.*F3 ( 12 )+Fr 
1  (  I?  )  ) /I ?.  +  3**( F3 ( II ) -3.*F3 ( 12 )+R.*F3 ( I  3 )-F2 ( 14) ) /R. ) 

F3 ( I )sF3 ( 11 )+E*( F4( II  )+0.5*( F4( II )-F4 (  12)  ) +5 • * ( F4 ( I]  ) -2 • * F4 ( 12  >+^4 
1(  13)  )/12«+3«*(F4( II ) — 3 • *F4 (  I  2 ) +3 • *F4 ( I  3 ) -F4 ( 1 4 ) ) /8 • ) 

F4  (I  )  =F4(  II  )  +F*  (  F5  (  I  1  )+0.5*(  FM  n  )-FF  (  I?  >  )+«S.*(  F*v(  1 1  )-2.*F5  (  12  >+F6 
1  (  13  )  )  /  I  2. +3.* ( F5 ( T 1  >-3.*F5 ( I  2  > +3.*F5 ( I  3 ) -F5 ( I  4 )  ) /R • ) 

F5 ( I ) =4 • *H*F  2 ( I ) *F 1 ( I )-(H+3. )*F4(  I ) *F0(  I  ) 

11  WRITE  (6*53)XI*F0(  I  )  *F1( I ) *F2( I ) *F3( I )  *F4( I )  *F5(  I ) 

IF (ABS(F3( I ) ) .LT.0.01 )  GO  TO  10 
1  =  1  +  1 
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MN=  I 

GO  TO  15 
10  CONTINUE 
STOP 

50  FORMAT ( 2X*23HTFMPFRATURF  EXPONENT  =  ,F8.4/) 

FORMAT(5X*2HXI*9X*2HFO*12X*2HF1>12X*?HF2*12X*2HF3*12X*2HF4*12X*2HF 

1  5  ) 

53  FORMAT ( 2X*F8»3*6F14*8 ) 

51  FORMAT  (  8F1.0. 6  > 

^ND 

SFNTRY  VFLTFM 
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APPENDIX  C 


DATA  REDUCTION 

The  typical  procedure  used  in  reducing  the  experi¬ 
mental  data  is  illustrated  by  considering  a  test  in  gly¬ 
cerine  for  w/6r  =  14.22.  Firstly,  all  thermocouple 
readings  were  converted  to  temperatures  using  the  conver¬ 
sion  chart  for  iron-constantan  thermocouples.  The  surface 
temperatures,  0  ,  were  then  plotted  versus  (x/w)  on  log- 
log  paper  and  the  temperature  index  n  was  determined  from 
the  slope  of  this  straight  line.  Since  in  most  cases 
there  was  not  a  thermocouple  located  exactly  at  the  root, 
0R  was  evaluated  by  extending  the  plot  back  to  x/w  =  1. 

For  the  example  under  consideration,  n  =  0.363  +  9.8%  and 
0R  =  12 . 3  5°F  +  1.2%. 

In  order  to  evaluate  the  theoretical  n  from  figure 
2.3  it  was  necessary  to  determine  the  numerical  value  of 
X»  For  evaluating  this,  a  bulk  fluid  temperature  was 
required  since  several  fluid  properties  were  needed.  The 


average  surface  temperature,  ®oavg*  Was  ta^en  as 


(--/  X)  ,  which  gave  0Qavg  =  9.06°F  +  11.0%. 


w 


This  along  with  T^  -  67.27°F  resulted  in  a  bulk  fluid  tempera¬ 
ture  of  Tbulk  =  7 1 . 80°F  +  11.5%. 

Properties  were  obtained  from  tables 


14 


,  giving 


. 


-r  x 
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a  =  11.65  x  103  +  3.0%,  =  88  +  4.0%  1/°F  ft3, 

v 

kf  =  .165  +  0.5%  BTU/hr . ft °F  and  k  =  132  +  2.0%  BTU/hr . f t . °F . 
Using  these  values, 


X 


req  9r  w  1 

\v 

-  /c2(l+a)  - 

1/4 


=  0.700  +  6.0%. 


Considering  the  experimental  temperature  profile  for 
a  typical  point  in  the  boundary  layer, 


x 

Y 

QU,y) 

w 


2"  +  0.05 
0.095”  +  0.002 
4 . 82°F  +  0.1 
6.938"  +  0.002 


6q  9R 

-4 k2  wn(l+a)- 


1/4 


_X_ 

1-n 


0.467  +6.1% 


<&(?) 


yR  W 


0.612  +5.7% 


These  percentage  errors,  though  calculated  for  the 
central  region  of  the  thermal  boundary  layer  in  glycerine, 
are  indicative  of  the  magnitude  of  the  errors  for  all 
tests  performed. 
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APPENDIX  D 

AN  ALTERNATIVE  SOLUTION 


During  the  course  of  attempting  to  solve  the  govern¬ 
ing  fifth  order  non-linear  equation  (9) ,  several  other 
numerical  and  analytical  methods  were  triedc  One  of  the 
numerical  methods  worth  noting  is  a  combination  of  numeri¬ 
cal  integration  and  iteration „  The  Fortran  IV  computer 
program  used  for  this  method  will  follow  the  outline  of 
the  procedure  used* 

The  uncombined  form  of  equation  (9) , 

F  ”  '  °  +  $  =  0  (Dl) 

0 '  '  +  (n  +  3)  F$  -  4n  f'$  =  0  (D2) 


is  the  starting  point 0  Formally  solving  (D2)  as  a  linear 

9 

non-homogeneous  first  order  equation  in  0  results  in 


0  (£) 


e 


(n+3)  Fd£ 


?  . 
4n  F  0 


0 


4 


(n+3)  Fd£ 


d£ 


(n+3) F  d£ 


’o 


+  $  (0)  e 


© 


(D3) 


Integrating  (D3)  once  and  using  the  boundary  condition, 


+  0  ( 0) 


•5 


(n+3)  Fd| 


e 


di  . (D4) 


1  ^  srii  :  II  i  j  b  s  e  o3  hs bu  :q 
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Applying  the  boundary  condition,  0(oo)  =  0,  results 


in 


oo 


oo 


4>  (0)  = 


p  oo 

-,1+  L  * 


(n+3)  Fd£ 


oo 


(n+3)  Fdf* 


4n  F  0  e 


oo 


(n+3)  Fd^ 


.00  uo 

e 


d£ 


at 


which  is  the  temperature  gradient  at  the  fin  surface . 

The  procedure  begins  by  assuming  an  approximate  temp¬ 
erature  profile,  $(§),  for  the  case  n  =  0.  Since  from 

•  •  •  ii 

(Dl) ,  F  -  -0,  the  corresponding  function  F  (£)  may  be 
found  by  numerically  integrating  $(£)  and  using  the 
appropriate  boundary  conditions  to  give 


I  9 


.co 


F  (?)  =  <t>(?)  d?. 

J? 

9 

The  functions  F  (§)  and  F(£)  are  obtained  in  the 

i  i  i 

same  manner  by  integrating  F  (£)  and  F  (£)  respectively 
Thus , 

9  P  ^  II 

F  (?)  =  /  F  (?)  d? 


F  (?) 


?  i 

F  (?)  d? 


(D6) 


(D7) 


Several  of  the  integrals  involved  require  integration 


(D5 


from  0  to  oo . 


It  is  obvious  from  figure  202  that  integration 


- 


*  ax  banxe^cfo  3Xb  {$)%  bae  ())  *3  znoilonui  srfT 


■ 
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from  £  =  0  to  3  is  sufficient  for  these  integrals.  All 
integrals  are  evaluated  using  the  trapezoidal  rule  and 
increments  of  0.05. 

Once  all  the  functions  and  integrating  factors  have 

0 

been  evaluated  $  (0)  may  be  determined  from  (D5)  and 
then  the  new  temperature  distribution ,  §(£),  evaluated 
from  (D4) .  It  is  now  necessary  to  check  this  new 
distribution  with  the  original  assumed  $(£)  for  each 
numerical  value  of  £9  which  goes  from  0  to  3  in  increments 
of  0.05.  Obviously  if  the  assumed  distribution  were 
correct  there  would  be  no  difference  between  the  two 
profiles.  However  if  for  any  value  of  £  considered,  the 
two  distributions  vary  by  more  than  the  set  limit  of 
0.001,  the  new  values  of  $(£)  are  now  used  as  the  start¬ 
ing  profile  and  the  procedure  is  repeated.  This  type 
of  iteration  is  carried  out  until  the  difference  between 
starting  and  final  profiles,  ,  for  every  value  of  £ 

considered,  is  less  than  0.001. 

When  used,  this  procedure  gave  temperature  and 
velocity  profiles,  for  values  of  n  up  to  0.4,  almost 
identical  to  those  presented  in  the  thesis.  For  larger 

values  of  n  difficulty  was  encountered  with  the  exponen- 

r? 

tial,  EXP  /  (n+3)Fd£,  since  the  exponent  became  larger 

'-'o 

than  that  allowed  in  the  computer.  However  it  is  quite 
possible  the  exponential  may  be  split  into  parts  and  this 
problem  overcome.  This  procedure,  though  not  complete,  is 
included  here  because  of  its  potential  as  a  method  for 
solving  problems  of  this  type. 
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SIBFTC  ALTSOL  NOL I  ST  ,NODFCK 

DIMENSION  F( 200  )  *F3 (3  *200) ,SUM( 200 )  *F2 ( 200 ) *Fl ( 200 ) *FO (20O ) ,0(20 ) 
1 *00( 200 ) *0R ( 200 ) »0S( 200 ) »0T ( ?00 ) *0U ( 200 ) *0V( 200 ) *01 ( 200) 

DIMENSION  U ( 200  )  ,UA ( 5  00 ) 

M=  1 0 

N  =  80 

TT=N 

EE=0.05 

F ( 1 )=1. 000000 

DO  1  J  =  2  ,  N 

P  =  J 

X= (P-1. ) /20. 

Y=X**1 • 4 

1  F ( J ) = .9530* EXP ( — Y ) 

DO  30  K= 1 *M 

R  =  K 

H= (  R  —  1  #  ) /10. 

WR ITE(6*50)H 
LL  =  0 
G=  ]  .0 
1=2 

11=1-1 
DO  2  J=1 * N 

2  F3 ( I  * J ) *F ( J ) 

20  DO  41  J= 1 *N 

41  F 3 (  1 1  * J ) =F3 ( I  * J ) 

C  F 3  IS  EQUIVALENT  TO  TEMP,  NOT  NFG  TEMP 

1  =  1-1 

SlJMM  =  0. 000000 
51JM  (  1  )  =0  •  000000 
DO  3  J  =  2  *  N 
J1=J-1 

SUM ( J ) = ( F3 ( I*J1 )+F3( I *J) )*EE 

3  SUMM=SUMM+SUM(  J > 

F2 ( 1 ) =SUMM 

C  F 2 ( 1 )  I S  THE  VELOCITY  GRADIENT  AT  THE  WALL 

DO  4  J=2*N 
SUMM=SUMM-SIJM(  J  ) 

4  F2 ( J ) =SUMM 
SUMM=0. 000000 
FI  ( 1  )=0,000000 
DO  5  J  =  2  ,  N 
J1=J-1 

SUM ( J ) = ( F2 ( J 1 )+F2( J) )*FE 
SUMM=SUMM+SUM ( J ) 

5  F 1 ( J ) =SUMM 
SUMM=0# 000000 
FO ( 1 ) =0.000000 
DO  6  J  =  2  *N 
J1=J-1 

SUM( J)=(F1( J1)+F1( J) ) *EE 
SUMM=SUMM+SUM( J ) 

6  FO(J)=SUMM 

C  Q(J)  IS  THE  INTEGRAL  OF  THF  INTEGRATING  FACTOR 

Q (  1  )  =0.000000 
FUMM=0. 000OO0 
DO  7  J=2*N 
J1=J-1 

SUM (J)=(H+3.)*(F0(J1>+F0(J) )*FE 
SUMM=SUMM+SUM ( J ) 
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3^“un+‘  fiu^^uMue 

(  L  J'OP.^MMUPxMmijp. 
VM(|?r|  L  ) 
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Pi 
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0( J ) -SUMM 
DO  1 R  J  =  1  ,N 

01  ( J ) =Q ( J ) / 1 0  • 

00(J)  IS  VFL*TEMP* INTEGRATING  FACTOR 
DO  8  J=1 »N 

OQ (J)=F3(I*J)*F1(J)*( {EXP(Q1(J) ) )**lo) 

OR ( J )  IS  THE  INTEGRAL  OF  QQ(J) 

OR  n ) =n . oooooo 
<;!  iMM=n,  oooooo 

DO  9  J  =  2  *N 
J1=J-1 

SUM ( J ) = ( OQ ( J 1 ) +  QQ ( j ) ) *FF 
SUMM  =  Sl)MM+SUM  (  J  ) 

OR ( J )=SUMM 

QS(J)  is  OR ( J ) *  INVFRSE  OF  INTEGRATING  FACTOR 
OS (l ) =0.000000 
DO  10  J  =  2  *N 

QS ( J ) =0R ( J ) * ( ( FXP ( -Ql ( J ) ) ) **i o ) 

OT ( J )  IS  THE  DOUBLE  INTEGRAL 
OT (1 ) =n#ooonoo 
SUMM=0. 000000 
DO  11  J  =  2  *N 
J1=J-1 

SUM( J)=4.*H*(QS( J1 )+QS( J) ) *EE 
SUMM=SUMM+SUM<  J ) 

QT( J)=SUMM 
DO  12  J=1 *N 

QU (  J  )  =  ( (FXP (-01 ( J ) ) ) ** l 0 ) 

SUMM=0. OOOOOO 
0V( 1 ) =0.000000 
DO  13  J  =  2  » N 
J1=J-1 

SUM ( J ) = ( OU ( J 1 ) +QU ( J ) ) *EE 
SUMM=SUMM+SUM( J ) 

QV( J)=SUMM 

OV(J)  IS  THE  INTEGRAL  OF  THE  INVERSE  OF  INT  FACTOR 
C*s—  (  1  .00  0000+0  T  (N)  )  /QV  (  N  ) 

C  IS  THE  TEMP  GRADIENT  AT  THE  WALL 
1=1  +  1 
11=1-1 
DO  16  J=1 >N 

F3 (  I  ,J  )  =1 .000000+  C*QV ( J ) +0T ( J ) 

IF(LL.GT.O)  GO  TO  AO 
L  =  1 

SUMM=0. 000000 
DO  21  J=  1  *N 

U( J)=0.01*( ABS(E3(  I  * J )-F3( I  1  * J)  )  ) 

SUMM=SUMM+U( J) 

UA ( L ) =SUMM/TT 
DO  22  J= 1 *N 

F3(I*J)=F3(I*J) +G*U ( J ) *F3 ( I *  J  ) 

L  =  L  + 1 
LL  =LL+ 1 
GO  TO  20 
SUMM=0. 000000 
DO  31  J  =  1  *  N 

U(J)=0.01*(ABS(E3 (  I , J )-E3( I  1 , J)  )  ) 

SUMM=SUMM+U( J> 

UA  (  L  )  srSUMM/TT 
Z=ABS(UA (L  )  ) 
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W=  ARS ! UA ( L-l )  ) 

IF (Z.GT.W>G=-G 

IF!7.LT. 0.000001 )  GO  TO  15 
DO  23  J= 1 *N 

2  3  F5  ( I  * J ) =F3 ( T * J )+G*U( J )*F3 ( I  *  J ) 

L  =  L+1 
GO  TO  20 
15  CONTINUF 

WR I TE ( 6  *  5 1 )C 
DO  14  J  =  1  *N 
S  =  J 

XI= (S-l . ) /20. 

14  WR  I TF ( 6 ♦ 53 )  X  I  *  FO ( J ) *  F 1 ( J ) *  F  2 ( J ) ♦ F  3 ( I ♦  J  ) 

DO  30  J=1*N 
30  F(J)=F3( I  *  J  ) 

FT  OP 

50  FORMAT! ?X*23HTEMPERATURE  EXPONENT  =  ,F8.4/) 

51  FORMAT! 2X*24HTEMP  GRADIFNT  AT  WALL  =  *F14.8/> 
53  FORMAT !2X*F8.4»4F20.8) 

55  FORMAT! 2X*8E16. 8 ) 

FND 
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